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ABSTRACT 

We investigate further a model of the accreting millisecond X-ray pulsars we proposed earlier. In 
this model, the X-ray-emitting regions of these pulsars are near their spin axes but move. This is 
to be expected if the magnetic poles of these stars are close to their spin axes, so that accreting gas 
is channeled there. As the accretion rate and the structure of the inner disk vary, gas is channeled 
along different field lines to different locations on the stellar surface, causing the X-ray-emitting areas 
to move. We show that this "nearly aligned moving spot model" can explain many properties of the 
accreting millisecond X-ray pulsars, including their generally low oscillation amplitudes and nearly 
sinusoidal waveforms; the variability of their pulse amplitudes, shapes, and phases; the correlations in 
this variability; and the similarity of the accretion- and nuclear-powered pulse shapes and phases in 
some. It may also explain why accretion-powered millisecond pulsars are difficult to detect, why some 
are intermittent, and why all detected so far are transients. This model can be tested by comparing 
with observations the waveform changes it predicts, including the changes with accretion rate. 
Subject headings: pulsars: general — stars: neutron — stars: rotation — X-rays: bursts — X-rays: 
stars 



, 1. INTRODUCTION 

Highly periodic millisecond X-ray oscillations have 
been detected with high confidence in 22 accreting neu- 

■ tron stars in low-mass X-ray binary systems (LMXBs), 
using th e Rossi X-ray Timing Exp lorer (RXTE) satel- 
lite (see lLamb &; Boutlou kos 2008|) . We refer to these 

OO ' stars as accreting millisecond X-ray pulsars (AMXPs). 
Accretion-powered millisecond oscillations have so far 
been detected in 10 AMXPs. They are always observ- 

■ able in seven AMXPs, but are only intermittently de- 
^ ' tected in three others. Nuclear-powered millisecond os- 

. t-H ! dilations have been detected with high confidence during 
thermonuclear X-ray bursts in 16 AMXPs. Persistent 
accretion-powered millisecond oscillations have been de- 
| tected in two AMXPs that produce nuclear-powered mil- 
lisecond oscillations; intermittent accretion-powered mil- 
lisecond oscillations have been detected in two others. 
The AMXPs have several important properties: 
Low oscillation amplitudes. The fractional am- 
plitudes of the accretion-powered oscillations of 
most AMXPs are often only -l%-2%. 1 Persis- 
tent accretion-powered oscillations with amplitudes 
< 1% are often det ected with hig h con fi dence in 
IGR J00291+5934 (IGallowav et all [20051: IPatrunol 
20081) and X TE J1751-305 (|Markwardt et alJ I2002T 
Patruno 2008) . Persistent accretion-powered oscillations 
with amplitudes as low as 2% are regularly seen in 

1 We characterize the strengths of oscillations by their rms ampli- 
tudes, because the rms amplitude can be defined for any waveform, 
is usually relatively stable, and is closely related to the power. We 
convert reported semi-amplitudes of purely sinusoidal oscillations 
or Fourier components to rms amplitudes by dividing by \/2. 



XTE J 1807- 294 (jZhang et al.1 [2fM IChou et al l [20081 : 
IPatrunol 12001 . XTE J092 9-314 (IGallowav et al.| [2002). 
and XTE J1814-338 (IChung et alJ 120081: IPatrunol 
2008). The amplitude of the accretion-powered os- 
cillation seen in SWIFT 1756.9-2508 was -6% 
(jKrimm et alJ 120071 ). The intermittent accretion- 
powered oscil l ations detected in SAX J 1748. 9-2021 
ilGavriil et all 120071 : lAltamirano et alJ [20081: IPatrunol 
[2008D, HETE J1900.1-245 5 (IGallowav et all I2007TT 
and Aql X-l (|Casella et alJ 120081 ) have amplitudes 
~0.5%-3%. 

Nearly sinusoidal waveforms. The waveforms (light 
curves) of the accretion-powere d oscillations o f most 
AMXPs are nearly sinusoidal (see Wiinands 2006 and the 
references in the preceding paragraph). The amplitude 
of the first harmonic (fundamental) component is usually 
> 10 times the amplitude of the second harmonic (first 
overtone) component, although the ratio can be as small 
as ~3.5, as is som etimes the case in XTE J1807— 294 
([Zhang et al.l l2006L or even > 1, as is sometimes the ca se 
in SAX J1808.4-3658 (see, e.g.. lHartman et al.|[2008l) . 

Highly variable oscillation amplitudes. The fractional 
amplitudes of the accretion-powered oscillations of most 
AMXPs vary in time by factors ranging from ~ 2 to 
~ 10. Observed fractional amplitudes vary from 0.7% 
to 1.7% in SAX J1748.9-2021, from 0.7% to 3.7% in 
XTE J1751-305, from 3% to 7% in XTE J0929-314, 
from 1% to 9% in IGR J00291+5934, from 2% to 11% in 
XTE J1814-338, from 1% to 14% in XTE J1808-338, 
and from 2% to 19% in XTE J1807-294 (see the refer- 
ences above). 

Highly variable pulse phases. The phases of accretion- 
powered pulses have been seen to vary rapidly by 
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as much as ~ 0.3 c ycles in several AMXPs, including 
SAX J1808.4-3658 tiMorgan et al.ll2003t lHartman et all 
[20081) and XTE J1807-294 (jMarkwardti I2004D . Wild 
changes in the apparent pulse frequency have been ob- 
served with both sig ns at the same a ccretion rate in 
XTE J1807-294 (see iMarkwardtl \2004) . If interpreted 
as caused by changes in the stellar spin rate, these phase 
variations would be more than a factor of 10 larger than 
expected for the largest accretion torques and smallest 
inertial moments thought possible for these system s (see 
IGhosh fc Lamblll979b| -" lLattimer fc Prakashll200lD . 

Undetected accretion-powered oscillations. More than 
80 accreting neutr o n stars in LM XBs are known 
(|Chakrabartvl l2005t iLiu et all l2007t ). but accretion- 
powered millisecond X-ray oscillations have so far been 
detected in only 10 of them. Accretion-powered oscil- 
lations have not yet been detected even in 13 AMXPs 
that produce periodic nuclear-powered millisecond oscil- 
latio ns, indicating that they ha ve millisecond spin peri- 
ods (|Lamb fc Boutloukosll2008D ; eight of these also pro- 
duce kilohertz quasi-periodic oscillations (QPOs) with 
frequency separations that indicate that they not only 
have millisecond spin perio ds but also have dynamica lly 
important magnetic fields (jBoutloukos fc Lamb 2008). 

Intermittent accretion-powered oscillations. Accretion- 
powered millisecond X-ray pulsations have been 
detected only occ a sionally in SAX J1748. 9-2021 
(IGavriil et all 120071 : lAltamirano et al.l l200l IPatrunol 
2008), HETE J1900.1-2455 (G alloway et al.ll2007D . and 
Aql X-l (jCasella et al.ll2008[) . When oscillations are not 
detected, the upper limits are typically < 0.5%. 

Correlated pulse arrival times and amplitudes. The 
phase residuals of the accretion-powered pulses of several 
AMXPs appear to be anti-correlated with their fractional 
amplitudes, at least over some of the amplitude ranges 
observed. AMXPs that show this type of beh avior in- 
clude XTE J1807-294 and XTE J1814-338 (jPatrunol 
12008ft . 

Similar accretion- and nuclear-powered pulses. The 
shapes and phases of the nuclear-po wered X-ray pulses 
of th e AMXPs SAX J1808. 4-3658 (IChakrabartv et all 
12001 and XTE J1814-338 (jStrohmaver et al.H2003D~ are 
very similar to the shapes and phases of their accretion- 
powered X-ray pulses. 

Concentration in transient systems of AMXPs with 
accretion-powered oscillations. The AMXPs in which 
accretion-powered oscillations have been detected tend 
to be found in binary sy stems that have outbu rsts lasting 
about a month (but see lGallowav et al. 2008) separated 
by quiescent inter vals lasting years (jChakrabartvl 120051 : 
Riggio et al. 2008). The accretion rates of these neutron 
stars are very low. 

In this paper we investigate further the "nearly aligned 
moving spot" model of AMXP X-ray emission that w e 
proposed previously (jLamb et all 120061 . 120071 l2008afl . 
This model has three main features: 

1. The strongest poles of the magnetic fields of neu- 
tron stars with millisecond spin periods are located 
near and sometimes very near the stellar spin 
axis. This behavior is expected for several mag- 
netic field evolution mechanisms. 

2. The star's magnetic field channels accreting gas 



close to its spin axis, creating X-ray emitting ar- 
eas there and depositing nuclear fuel there. 2 

3. The X-ray emitting areas on the stellar surface 
move, as changes in the accretion rate and the 
structure of the inner disk cause accreting gas to be 
channeled along different field lines to different lo- 
cations on the stellar surface. (The magnetic field 
of the neutron star is fixed in the stellar crust on 
the timescales relevant to the phenomena consid- 
ered here.) 

These features provide the basic ingredients needed to 
understand the AMXP properties discussed above. This 
is the subject of the sections that follow. As a guide to 
these sections, we summarize our results here. 

1. Emission from near the spin axis naturally pro- 
duces weak modulation, regardless of the viewing 
direction. The reason is that uniform emission from 
a spot centered on the spin axis is axisymmctric 
about the spin axis and therefore produces no mod- 
ulation. Emission from a spot close to the spin axis 
has only a small asymmetry and therefore produces 
only weak modulation. 

2. Emission from near the spin axis also naturally pro- 
duces a nearly sinusoidal waveform, because the 
asymmetry of the emission is weak and broad. 

3. If the emitting area is close to the spin axis, even 
a small movement in latitude can change the oscil- 
lation amplitude by a substantial factor. 

4. If the emitting area is close to the spin axis, move- 
ment in the longitudinal direction by a small dis- 
tance can change the phase of the oscillation by a 
large amount. 

Changes in the latitude and longitude of the emit- 
ting area are expected on timescales at least as 
short as the ^0.1 ms dynamical time at the stellar 
surface and as long as the ~ 10 d timescale of the 
variations observed in the mass accretion rate. 

5. If the emitting area is very close to the spin axis 
and remains there, the oscillation amplitude may 
be so low that it is undetectable. The effects of 
rapid changes in the position of the emitting area — 
possibly in combination with other effects, such 
as reduction of the modulation fraction by scat- 
tering in circumstellar gas — may also play a role 
in reducing the detectability of accretion-powered 
oscillations in neutron stars with millisecond spin 
periods. These effects may explain the fact that 
accretion-powered X-ray oscillations have not yet 
been detected in many accreting neutron stars that 
are thought to have millisecond spin periods and 
dynamically important magnetic fields. 



2 IMuno et a l. ( 2002) considered a single bright spot near the spin 
axis as well as a uniformly bright hemisphere and antipodal spots 
near the spin equator as possible reasons for the nearly sinusoidal 
waveforms of some X-ray burst oscillations, but did not consider 
accretion-powered oscillations or other consequences of emission 
from near the spin axis. 
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6. If the emitting area is very close to the spin axis, 
a small change in the accretion flow can suddenly 
channel gas farther from the spin axis, causing the 
emitting area to move away from the axis. This can 
make a previously undetectable oscillation become 
detectable. Temporary motion of the emitting area 
away from the spin axis may explain the intermit- 
tent accretion-pow ered oscillations of some AMXPs 
(jLamb et al.|[2009h . 

7. If the pulse amplitude and phase variations ob- 
served in AMXPs are caused by motion of the emit- 
ting area, they should be correlated. In partic- 
ular, the pulse phase should be much more scat- 
tered when the pulse amplitude is very low. The 
reason is that changes in the longitudinal position 
of the emitting area by a given distance produce 
much larger phase changes when the emitting area 
is very close to the spin axis, which is also when 
the oscillation amplitude is very low. 

The observational consequences discussed so far de- 
pend only on features (2) and (3) of the model, 
i.e., that the accretion-powered X-ray emission of 
AMXPs comes from areas near their spin axes and 
that these areas move significantly on timescales of 
hours to days. 

8. The picture of AMXP X-ray emission outlined here 
suggests that the shapes and phases of the nuclear- 
and accretion-powered pulses are similar to one an- 
other in some AMXPs because the nuclear- and 
accretion powered X-ray emission comes from ap- 
proximately the same area on the stellar surface. 
The reason for this is that in some cases, the mecha- 
nism that concentrates the magnetic flux of the ac- 
creting neutron star near its spin axis, as it is spun 
up, will naturally produce magnetic fields strong 
enough to confine accreting nuclear fuel near the 
magnetic poles at least partially, even though the 
dipole component of the magnetic field is weak. 

9. The picture of neutron star magnetic field evolution 
and AMXP X-ray emission outlined here also sug- 
gests a possible explanation for why the AMXPs 
in which accretion-powered oscillations have been 
detected arc in transient systems. If most neu- 
tron stars in LMXBs were spun up by accretion 
from a low spin rate to a high spin rate, their mag- 
netic poles were forced very close to their spin axes, 
making accretion-powered oscillations difficult or 
impossible to detect. However, those stars that 
are now in compact transient systems now expe- 
rience infrequent episodes of mass accretion and 
the accretion rate is very low. By now they have 
been spun down from their maximum spin rates, a 
process that could force their magnetic poles away 
from their spin axes enough to produce detectable 
accretion-powered oscillations. 

These last two observational consequences depend on 
feature (1) of the model, i.e., on how the magnetic fields 
of neutron stars evolve as they are spun up and down by 
accretion and electromagnetic torques. 

In the remainder of this paper we discuss in detail the 
features of the model and its observational implications. 



In Section [21 we outline our approach, discussing our 
modeling of X-ray emission from the stellar surface, our 
computational and the code verification methods, and 
the pulse profile representation we use. 

We present our results in Sections [3] and [U These re- 
sults are based on our computations of several hundred 
million waveforms for different emitting regions, beam- 
ing patterns, stellar models, and viewing directions. In 
Section [3l we consider the shape and amplitude of X- 
ray pulses as a function of the size and inclination of the 
emitting areas, the compactness of the star, and the stel- 
lar spin rate. In SectionlH we consider the changes in the 
pulse amplitude and phase produced by various motions 
of the emitting regions on the stellar surface and explore 
the origins of correlated changes in the pulse amplitude 
and phase and the effects of rapid movement of the emit- 
ting areas. We also discuss why oscillations have not yet 
been detected in many accreting neutron stars in LMXBs 
and why accretion-powered oscillations are detected only 
intermittently in some AMXPs. 

In Section [5l we summarize the results of our model 
calculations. We also discuss how the magnetic poles of 
most AMXPs can be forced close to their spin axes, how 
such mechanisms may explain why the AMXPs that pro- 
duce accretion-powered millisecond oscillations are tran- 
sient pulsars, the consistency of the model with the ob- 
served properties of rotation-powered millisecond pul- 
sars, and possible observational tests of the model dis- 
cussed here. 

Further results of our investigation of the present 
model will be presented elsewhere (S. Boutloukos et al., 
in preparation). 

2. X-RAY WAVEFORM MODELING 

2.1. Modeling the X-ray Emission 

In the radiating spot model of AMXP X-ray emission, 
the waveform seen by a distant observer depends on the 
sizes, shapes, and positions of the emitting regions on the 
stellar surface; the beaming pattern of the radiation; the 
compactness, radius, and spin rate of the star; and the 
direction from which the star is observed. The proper- 
ties of the X-ray emitting regions are determined by the 
strength and geometry of the star's magnetic field, the 
locations where plasma from the accretion disk enters the 
magnetosphere, the extent to which the accreting plasma 
becomes threaded and channeled by the stellar magnetic 
field, and the resulting plasma flow pattern onto the stel- 
lar surface. 

In principle, accreting plasma can reach the stellar 
surface in two basic ways: (1) by becoming threaded 
by the stellar magnetic field and then guided along 
stellar field lines to the vic i nity of a stellar mag- 
netic pole (lLamb et al.1 119731: iBasko fc Sunvaevl 119751: 



Eisner fe Lainbl lT976: Ghosh et al. 1 119771 : IGhosh fc Lambl 
1979allbD or (2) by penetrating between lines of the stellar 
magnetic field via the magnetic ve rsion of the Rayleigh- 
Taylor instability (|Lambljl975allb1: lElsner fc Lamblll976l 
1977t lArons fc Lealll976t lLamblll977h and then spiraling 
inward to the stellar surface. 

If a centered dipole component is the strongest compo- 
nent of the star's magnetic field, plasma in the accretion 
disk that becomes threaded and then channeled to the 
vicinity of a magnetic pole is expected to impact the star 
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in a partial or complete annulus around the pole, produc- 
ing a crescent- or ring-shaped emitting area near the pole. 
If the axis of the dipole field is significantly tilted rela- 
tive to the spin axis and the spin axis is aligned with the 
axis of the accretion disk, a crescent-shaped emitt i ng re- 
gion is expected (s e e. e.g.JBasko fc Sunvaeylll97. ' 



Ghosh et al 



Miller et al.l 



119771: iDaumerie et al.1 



Sunvacv 1975. 
119961: IMilleif 



1976 



1996 



19981 : IRomanova et al.1 12003) . If instead the simultaneously (see 



l200l . 

Whether accreting plasma reaches the stellar surface 
predominantly via channeled flow along field lines or via 
unstable flow between field lines depends on the accre- 
tion r ate and the spin frequency of the star (see iLambl 
[1981 iRomanova et all [20081 : IKulkarni et al.ll200l . Un- 
der some conditions , plasma may accrete in both ways 



dipole axis is very close to the spin axis, as in the model 
of AMXP X-ray emission proposed here, the emitting 
region may completely encircle the spin axis ( see, e.g., 
IGhosh fc Lamhlll979aU bl: IRomanova et^II[200l . 

The north and south magnetic poles of some AMXPs 
may be very close to the same spin pole, producing a 
ver y off-center dipole moment orthogonal to the spin axis 
(see lChen fe Rudermanll 19931: 1 Chen et al.lll998t and Sec- 
tion [5TT]) . If so, accreting matter will be channeled close 
to the spin axis, but may be channeled preferentially to- 
ward one magnetic pole, producing an emitting region 
with approximately one-fold symmetry about the spin 
axis, or about equally toward both poles, producing an 
emitting region with approximately two-fold symmetry 
about the spin axis. In the first case the first harmonic 
of the spin frequency is likely be the dominant harmonic 
in the X-ray waveform whereas in the second case the 
second harmonic is likely to dominate. Which case oc- 
curs will depend on the accretion flow through the inner 
disk. In either case, the X-ray emission will come from 
close to the spin axis. 

The neutron stars that are AMXPs may well have 
even more complicated magnetic fields, with significant 
quadrupole and octopole components. Higher multipolc 
components are likely to play a more important role in 
the AMXPs than in the classic strong-field accretion- 
powered pulsars, because the magnetic fields of the 
AMXPs are much weaker. As a result, accreting plasma 
can penetrate closer to the stellar surface, where the 
higher multipole moments of the star's magnetic field 
have a greater influence on t he channeling of accreting 
plasma ([Eisner fc Lamblll976[ ). In this case, plasma will 
still tend to be channeled toward regions on the surface 
where the magnetic field is strongest and will tend to 
impact the surface in rings or annuli, but the emission 
patte rn may be spatia lly complex and vary rapidly in 
time ([Long et al.ll2008ft . 

Disk plasma that penetrates between lines of the stel- 
lar magnetic field will continue to drift inward as it loses 
its angular momentum, probably predo minantly via its 
inter action with the star's magnetic field (|Lamb fc Miller 
2001). Cold plasma will remain in the disk plane 
and impact the star in an annul us where the disk 
plane intersects the stellar surface ( Miller et al.l 119981 : 
lLamb fc Miflerl [200lT ) . If some of the accreting plasma 
were to become hot, the forces exerted on it by the stel- 
lar magnetic field would tend to drive it toward the star's 
magnetic equator (Michel 1973), causing it to impact the 
stellar surface in an annulus around the star's magnetic 
equator. However, emission and inverse Compton scat- 
tcring of radiation is likely to keep the accreting plasma 
cold ([Eisner fc Lambl[l984f ). so that it remains in the disk 
plane as it drifts inward. Plasma that penetrates to the 
stellar surface via the magnetic Rayleigh-Taylor instabil- 
ity is likely to impact the stellar surface in rapidly fluc- 
tuating, irregular patterns (see IRomanova et al.1 120061 . 



2008; Kul karni et al 



a y accr ete 
Miller et all [19981 : [R 
120081) . 
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The sizes, shapes, and locations of the emitting ar- 
eas on the surface of an accreting magnetic neutron 
star and the properties of the emission from these 
areas are expected to change on timescales at least 
as short as the ~ 1 ms dynamical timescale near the 
star. This expectation is supported by recent simula- 
t ion s of accretion onto weakly magnetic neutron stars 
(seelRomanova et al.ll2003ll2004 12006[ lLong et al.l l2008t 
IRomanova et al.l 120081 : IKulkarni et al.l l2008h . However, 
changes in AMXP X-ray fluxes can be measured ac- 
curately using current instruments only by combining 
~ 100-1,000 s of data and hence only variations in wave- 
forms on timescales longer than this can be measured 
directly. Consequently, the emitting areas and beaming 
patterns that are relevant for comparisons with current 
observations of waveforms are the averages of the actual 
areas and beaming patterns over these relatively long 
times. The emitting areas and beaming patterns that we 
use in our computations should therefore be interpreted 
as averages of the actual areas and beaming patterns over 
these times. 

We have computed the X-ray waveforms produced by 
emitting regions with various sizes, shapes, and posi- 
tions, for several different X-ray beaming patterns and 
a range of stellar masses, compactnesses, and spin rates. 
We find that in many cases these waveforms can be ap- 
proximated by the waveforms generated by a circular, 
uniformly emitting spot located at the centroid of the 
emitting region. The main reason for this is that an ob- 
server sees half the star's surface at a time (or more, 
when gravitational light deflection is included), which 
diminishes the influence of the size and detailed shape 
of the emitting region on the waveform. Consequently, 
we focus here on the waveforms produced by uniformly 
emitting circular spots. We will discuss the waveforms 
produced by emitting areas with other shapes, such as 
rings or crescents, in a subsequent paper (S. Boutloukos 
et al., in preparation). 

In addition to studying the X-ray waveforms produced 
by emitting areas with fixed sizes, shapes, positions, and 
radiation-beaming patterns, we are also interested in the 
changes in waveforms produced by changes in the these 
properties of the emitting areas. The changes we investi- 
gate should be understood as occurring on the timescales 
> 100 s that can be studied using current instruments. 
It is not yet possible to compute from first principles the 
accretion flows and X-ray emission of AMXPs on these 
timescales, so simplified models must be used. (The sim- 
ulations referred to earlier follow the accretion flow for 
a few dozen spin periods or dynamical times, intervals 
that are orders of magnitude shorter than the intervals 
that are relevant). 

In the following sections, we consider radiation from a 
single spot, from two antipodal spots, and from two spots 
in the same rotational hemisphere, near the star's spin 
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axis. Although many uncertainties remain, recent mag- 
netohydrodynamic simulations of accretion onto weakly 
magnetic neutron stars ha ve found that gas impa cts \%- 
20% of the stellar surface (jRomanova et al.ll2004l ). equiv- 
alent to the areas of circular spots with angular radii of 
10°-53°. These radii are consistent with analytical esti- 
mates of the sizes of the emission regio ns of accreting neu- 
tron stars with weak magne tic fields (|Miller et al.lll998t 
iPsaltis fc Chakrabartvi[l999l ) . Consequently, we focus on 
spot sizes in this range. 

An observer may see radiation from a single spot either 
because the accretion flow pattern strongly favors one 
pole of a dipolar stellar magnetic field over the other, or 
because the observer's view of one pole is blocked by the 
inner disk or by accreting plasma in the star's magne- 
tosph ere (see iMcCrav fc Lambl 119761 : iBasko k, Sunyaevl 
1976). An observer may see radiation from two antipo- 
dal spots if emission from both magnetic poles is visible. 
Finally, an observer may see radiation from two spots 
near the same rotation pole if neutron vortex motion 
drives both of the star's di polar magnetic p oles toward 
the same rotation pole (see[Chenet23[r998). 

To make it easier for the reader to compare cases, we 
usually report results for our "reference" star, which is a 
1.4 Mq star with a radius of 5M in units where G = c = 1 
(10.3 km for M = 1.4 M©), spinning at 400 Hz as mea- 
sured at infinity, but we also discuss other stellar models. 
For the same reason, we usually consider spots with an- 
gular radii of 25°. This is not an important limitation, 
because the observed waveform depends only weakly on 
the size of the emitting spots, as discussed in Section l3~2l 
We describe how the results change if the spot is larger 
or smaller. 

2.2. Computing X-ray Waveforms 

The X-ray waveforms calculated here assume that radi- 
ation propagating from emitting areas on the stellar sur- 
face reaches the observer without interacting with any in- 
tervening matter. The bolometric X-ray waveforms that 
would be seen by a distant observer were calculated using 
the Schwarzs child plus Dopp l er (S+ D) approximation in- 
troduced bv iMiller & Lambl (|1998f ). The S+D approxi- 
mation treats exactly the special relativistic Doppler ef- 
fects (such as aberrations and energy shifts) associated 
with the rotational motion of the stellar surface, but 
treats the star as spherical and uses the Schwarzschild 
spacetime to compute the general relativistic redshift, 
trace the propagation of light from the stellar surface to 
the observer, and calculate light travel-time effects. For 
the stars considered here, and indeed for any stars that 
do not both rotate rapidly and have very low compact- 
ness, the effects of stellar oblateness and frame dragging 
are minimal and are neglig ible compared to un certainties 
in the X-ray emission (see iCadeau et al~ll2007ft . 

We describe the emission from the stellar surface us- 
ing coordinates centered on the star. When considering 
emission from a single spot, we denote the angle between 
its centroid and the star's spin axis by i s and its azimuth 
in the stellar coordinate system by <j) s . When consid- 
ering emission from two spots, we somewhat arbitrarily 
identify one as the primary spot and the other as the sec- 
ondary spot. We denote the inclination and azimuth of 
the centroid of the primary spot by i s \ and <j) s \ and the 
inclination and azimuth of the centroid of the secondary 



spot by i S 2 and <p 8 2- We denote the inclination of the 
observer relative to the stellar spin axis by i. 

In computing the waveforms seen by distant observers, 
we use as our global coordinate system Schwarzschild 
coordinates (r, 8, ip, t) centered on the star with 9 = 
aligned with the star's spin axis and ip = in the plane 
containing the spin axis and the observer. We choose 
the zero of the Schwarzschild time coordinate t so that a 
light pulse that propagates radially from a point on the 
stellar surface immediately below the observer (i.e., at 
8 = i and (p = 0) arrives at the observer at t = 0. 

We carried out many calculations to test and verify the 
computer code used to obtain the results we report here. 
We determined that the code was giving sufficiently accu- 
rate results by varying the spatial and angular resolutions 
used. For most of the cases considered in this paper, the 
emitting spots were sampled by a grid of 250 points in lat- 
itude and 250 points in longitude, the radiation-beaming 
pattern was specified at 10 4 angles, and the flux seen by 
a distant observer was computed at 10 4 equally spaced 
values of the star's rotational phase. In some cases, finer 
grids were used. 

We verified the code used here by comparing its re- 
sults with analytical and numerical results for several 
test cases: 

1. We tested our code's representation of emitting ar- 
eas and ray tracing in flat space by comparing the 
results given by our code with exact analytical re- 
sults for the absolute flux seen by an observer di- 
rectly above uniform, isotropically emitting circu- 
lar spots of various sizes. The numerical results 
agreed with the analytical results. 

2. We tested our code's computation of special rel- 
ativistic Doppler boosts, aberrations, and propa- 
gation time effects in several ways. We compared 
the results given by our code with exact analytical 
results for the waveforms produced by emission in 

(a) a pencil beam normal to the surface and (b) a 
thin fan beam tangent to the surface of a rapidly ro- 
tating star. We also compared the results given by 
our code with analytical results for the waveforms 
produced by a small spot on the surface of a slowly 
rotating star in flat space emitting (a) isotropically 
and (b) in a beamin g pattern representing Comp - 
tonized emission (see lPoutanen fc Gierlih ski 2003). 
The numerical results agreed with the analytical 
results. 

3. We tested our code's computation of the general 
relativistic redshift and light deflection for nonro- 
tating stars by (a) comparing the deflection of a 
fan beam tangent to the stellar surface given by 
our code for a variety of stellar compactnesses with 
the anal ytical expression s for t he li ght deflectio n 
given bv iPechenick et all (fl98l and [Pail (fl995T) : 

(b) comparing the absolute flux given by our code 
for an observer directly above isotropically emit- 
ting uniform circular spots of various sizes with in- 
dependent semi-analytical results for these cases; 

(c) comparing the symmetries of the waveform and 
the phase of the waveform maximum given by our 
code with exact analytical results for these quan- 
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tities; and (d) comparing the shape of the wave- 
for ms given by our code with the shapes r e porte d 
bv lPechenicket all (|1983f) and iStrohmaveil (fl992f) . 
Our numerical results agreed satisfactorily with the 
comparison results in all cases. 

4. We tested our code's computation of the wave- 
forms produced by emission from slowly rotating 
stars in general relativity by comparing the rms 
oscillation amplitudes it gives with the amplitudes 
given by the approximate a nalytical formulae of 
iViironen fe Poutanenl (|2004fl . Where the results 
of Viironen & Poutanen are expected to be ac- 
curate, the two sets of rms amplitudes agreed to 
better than 1%; in many cases the agreement was 
much better. We also compared the waveform 
given by our code for an isotropically emitting 
spot inclined 45° from the spin axis of a 1.4M Q 
star with a Schwarzschild coordinate radius of 5M 
spinning at 600 Hz seen by an observer at an in- 
clination of 45° with the wave form reported by 
iPoutanen fc Beloborodovl (2006) for this case; the 
two waveforms agreed to better than 1%. 

Further details of these tests and comparisons will be 
given in a subsequent paper (S. Boutloukos et al., in 
preparation). 

2.3. Constructing Pulse Profiles 

The X-ray flux seen by a given observer will evolve 
continuously in time as the star rotates and the emission 
from the stellar surface changes, generating the observed 
waveform W(t). As noted in Section |2~TI the accretion 
flow from the inner disk to the stellar surface is expected 
to vary on timescales at least as short as the ~ 1 ms 
dynamical timescale near the neutron star, which will 
cause the sizes, shapes, and positions of the emitting 
regions, and therefore the observed waveform, to vary on 
these timescales. 

The sensitivity of current instruments is too low to 
measure the waveform of an AMXP on timescales as 
short as 1 ms. However, nearly periodic waveforms with 
periods this short can be partially characterized by fold- 
ing segments of flu x data centered at a sequence of cloc k 
times U (see, e.g., IHartman etHl I2008t IPatrunol feOOS). 
If the data are folded with a period Pt that is chosen to 
agree as closely as possible with the local, approximate 
repetition period P(U) of the waveform, one can con- 
struct a time sequence of pulse profiles Wp(4>, ti); here <j> 
is the pulse phase over one cycle. 

The pulse profiles Wp{4>, U) constructed by folding flux 
data are averages of the actual pulse profiles over the time 
interval required to construct a stable profile, which can 

3 The waveforms reported by|Pcchcnick et al. (1983) for two an- 
tipodal spots are slightly inaccurate, as shown by the following 
two tests. The waveform seen by an observer in the star's rotation 
equator viewing two identical antipodal spots in the rotation equa- 
tor should be the same at 180° as at 0° , but this is not the case for 
their waveform for this geometry (see their Figure 7). More gener- 
ally, the flux from a uniform, isotropically emitting, circular spot 
on a spherical, nonrotating star should depend only on the angle 
between the radius through the center of the spot and the radius 
to the observer. This is n ot quite true for the waveforms reported 
by IPec hcnick et al. (1983). The waveforms given by our code pass 
these tests (for details, see S. Boutloukos et al., in preparation). 



be hundreds or even thousands of seconds, 10 5 -10 6 times 
longer than the ~ 1 ms dynamical timescale near the neu- 
tron star. The folded pulse profiles are therefore likely 
to vary more slowly and have less detail than the X-ray 
waveform, a point to which we will return in Section [3J 
The waveforms of AMXPs can be modeled even on the 
dynamical timescale near the stellar surface, but such 
waveforms would contain much more information than 
can be studied using current observations. Consequently, 
we focus here on modeling folded pulse profiles. We de- 
fine a computed pulse profile as the waveform seen by 
a given observer when a star with a constant emission 
pattern makes one rotation. It is customary and useful 
to describe pulse profiles by the amplitudes and phases 
of their Fourier components. We describe our computed 
profiles using the representation 

oo 

W P {t) =A + J2 A k cos[27rfc(i/ s i - <j) k )] . (1) 

k=l 

Here t is the time measured at infinity, A$ is the average 
of the flux over one rotation period, A)~ and are the 
amplitude and phase of the kth harmonic component of 
the pulse, and v s is the stellar rotation frequency. For 
harmonic k, the range of unique phases is to 1/k. We 
sometimes also refer to the phase <fi p of the peak of the 
pulse profile, defined in a similar way. With these defi- 
nitions, shifting the arrival time of a given pulse by an 
amount St shifts the phases of its Fourier components 
and the phase of the pulse peak by v s 8t cycles. 

We define the arrival time tk of the A:th harmonic com- 
ponent of the pulse profile to be zero if the first maximum 
of the sinusoid associated with this component reaches 
the observer at the same time as would a light pulse 
coming from a point on the stellar surface immediately 
below the observer (see Section 12. 2\ . The arrival time 
tk is affected by the sizes, shapes, and locations of the 
emitting areas, the beaming pattern of the emitted ra- 
diation, and the aberration and Doppler shifts produced 
by the spin of the star and usually is not zero. We de- 
fine the arrival phase of the kth harmonic component by 
4>k = v a tk ■ The arrival time t p and arrival phase 4> p of 
the peak of the pulse profile are defined similarly. We 
refer to tk and t p as time residuals and to 4>k and <p p as 
phase residuals. An increase in a time or phase residual 
indicates that the harmonic component or the peak is 
arriving later. 

3. OSCILLATION AMPLITUDES 

As discussed in Section [TJ the fractional amplitudes 
of the accretion-powered oscillations of most AMXPs 
are typically ~ l%-2%, but the amplitudes of several 
AMXPs are sometimes as large as ~ 10%-20%. A suc- 
cessful model of the accretion-powered oscillations of the 
AMXPs should therefore be able to explain oscillation 
amplitudes as low as ~ l%-2% without requiring a spe- 
cial viewing angle or stellar structure and should also be 
able to explain the higher amplitudes sometimes seen. 

As was also discussed in Section [JJ the accretion- 
powered oscillations of the AMXPs are often nearly si- 
nusoidal. The second harmonic of the fundamental os- 
cillation frequency has been detected in seven of the ten 
known AMXPs, but it is typically > 10 times weaker than 
the fundamental, although in a few cases it is not this 
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weak and in one case, SAX J1808.4— 3658, it is sometimes 
stronger than the fundamental. Upper limits ~0.1%- 
0.2% have been placed on the amplitude of any second 
harmonic in three AMXPs. 

Upper limits ~ 0.5%-l% have been placed on the frac- 
tional amplitudes of accretion-powered millisecond oscil- 
lations in several AMXPs that produce nuclear-powered 
millisecond oscillations, for frequencies close to the fre- 
quencies of the nuclear-powered oscillations. A model 
that can explain these nondetections as well as the ob- 
served properties of the detected accretion-powered os- 
cillations would be very attractive. 

In this section, we show that emitting regions on or 
near the stellar surface can produce oscillation ampli- 
tudes as low as ~ l%-2% for a substantial range of view- 
ing directions only if they are located within a few de- 
grees of the stellar spin axis. Regions near the spin axis 
also naturally produce nearly sinusoidal pulse profiles. If 
they are very close to the spin axis and remain there, 
the amplitudes of the oscillations they produce can be 
~ 0.5% or less. Thus, emission from very close to the 
spin axis, in combination with other effects, such as rapid 
phase and amplitude fluctuations and suppression of X- 
ray mo dulation by scattering of photon s in circumstellar 
gas (see lLamb et alj|1985t [Miller 200]J , may explain the 
nondetection of accretion-powered oscillations in accret- 
ing neutron stars in which nuclear-powered oscillations 
have been detected. 

We show further that although the fractional oscilla- 
tion amplitudes produced by large regions tend to be 
smaller, this is a weak effect. Unless almost the entire 
surface of the star is uniformly emitting, amplitudes as 
small as those observed are possible even for large spots 
only if they are located near the stellar spin axis. We 
also show that although the fractional amplitudes pro- 
duced by very compact neutron stars tend to be smaller 
than the amplitudes produced by less compact stars, 
this effect is too weak to explain the small amplitudes 
of the AMXPs. Furthermore, it cannot explain why 
the fractional amplitudes of several AMXPs are ~ 1%- 
2% at some times but ~ 15%-25% a few hours or days 
later, since the stellar compactness cannot change on 
such short timescales. Finally, we note that the fractional 
amplitudes produced by spinning neutron stars depend 
only weakly on their spin rates. 

Measured oscillation amplitudes are likely to be smaller 
than those shown in the figures in this section, which are 
for stable spots fixed on the stellar surface. The reason 
is that the pulse shape is expected to vary on timescales 
shorter than the time needed to construct a pulse wave- 
form. Such rapid pulse shape variations will appear as 
increased background noise, red ucing the apparen t am- 
plitude of the oscillations (see lLamb et al.l [1985) • As 
discussed below, this effect tends to be more important 
when the emitting area is near the spin axis, and it will 
therefore tend to reduce further the apparent amplitude 
of the oscillations produced by spots at small inclina- 
tions. 

In Section 2) we show that the variations of the oscil- 
lation amplitudes seen in all AMXPs and the larger frac- 
tional amplitudes ~ 15%-20% occasionally seen in sev- 
eral of them can be explained by modest increases in the 
inclinations of their emitting regions, if these regions are 



close to the star's spin axis. 

3.1. Dependence on Spot Inclination 

The precise distance emitting regions can be from 
the spin axis and still produce oscillation amplitudes as 
low as the ~ l%-2% amplitudes often observed in the 
AMXPs depends on the beaming pattern of the emis- 
sion. The angular pattern of the radiation flux produced 
by Thomson scattering is strongly peaked in the direc- 
tion normal to the stellar surface, whereas the flux pat- 
terns produc ed by fan-shaped radiati o n bea ms like those 
proposed by IPoutanen fc Gierlihskil (|2003f) to describe 
Comptonized emission from near the stellar surface are 
peaked much less (if at all) in the direction normal to 
the stellar surface. The angular pattern of the radia- 
tion flux produced by isotropic emission is intermediate 
between these two cases, being peaked normal to the sur- 
face less than the pattern produced by Thomson scatter- 
ing but more than the pattern produced by fan-beam 
emission. We first discuss the waveforms produced by 
isotropic emission and then the waveforms produced by 
other beaming patterns. 

Figure [T] shows the fractional rms amplitudes of the 
first harmonic (fundamental) and second harmonic (first 
overtone) components of the pulse profile and the total 
fractional rms amplitude produced by isotropic emission 
from a single stable spot and from two stable antipodal 
spots, as a function of the observer's inclination and the 
inclination of the primary spot relative to the spin axis. 
The curves are for spots with angular radii of 25° on our 
reference star (a 1.4M star with a radius of 5M spinning 
at 400 Hz). As expected, all amplitudes are zero when 
the star is viewed along its spin axis (i = 0°) or the emit- 
ting spots are centered on the rotation pole (i s = 0° or 
i s i = 0° and i S 2 — 180°). In addition, the amplitude of 
the first harmonic variation produced by two stable an- 
tipodal spots vanishes if either the observer or the spots 
are in the rotation equator (i s = 90° or i s \ — i S 2 = 90°), 
whereas the amplitude of the second harmonic is highest 
for these geometries. 

The top panels in Figure [1] show that if the observer 
sees radiation only from a single spot, the total frac- 
tional amplitude can be as low as ~ l%—2% for a range 
of viewing angles only if the spot is within a few degrees 
of the spin axis. These panels also show that for radi- 
ation from a single spot, the second harmonic is ~5— 10 
times smaller that the first harmonic for most observer 
and spot inclinations. 

The bottom panels in Figure Q] show that if the ob- 
server sees radiation from two antipodal spots, the total 
fractional amplitude can be as low as ~ 2% for a substan- 
tial range of viewing angles only if the spots are within 
about 5° of the spin axis. These panels also show that 
for radiation from two antipodal spots, the second har- 
monic is smaller than the first harmonic for observer or 
spot inclinations less than about 60°. The amplitude of 
the second harmonic is < 2% for all observers only if the 
spots are within 15° of the spin axis. 

Two emitting spots centered near the same rotation 
pole are expected if the inward motion of neutron vor- 
tices in the stellar core drives both of the star's magnetic 
poles toward the same rotation pole (see Section 15. Q . 
The fractional modulation produced by two spots with 
this geometry is typically somewhat larger than the mod- 
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Fig. 1. — Top panels (left to right): contour plots of the fractional rms amplitudes of the first harmonic (fundamental) and second 
harmonic (first overtone) components of the pulse profiles and the total fractional rms amplitude produced by a single stable spot, as 
a function of the observer's inclination and the inclination of the primary spot relative to the spin axis. Bottom panels (left to right): 
corresponding contour plots for the pulse profiles produced by two stable antipodal spots. All plots assume spot radii of 25° and a XAMq 
star with a radius of 5M spinning at 400 Hz. The heavier lines highlight the contours where the fractional flux variation is 2%. This figure 
shows that emitting regions on or near the stellar surface can produce oscillation amplitudes as low as ~ 1%— 2% for a substantial range of 
viewing directions only if they are located within a few degrees of the stellar spin axis. 



ulation produced by two antipodal spots if the observer 
and spot inclinations are both large, but is typically 
smaller than the modulation produced by two antipo- 
dal spots if the observer inclination or the spot inclina- 
tion is small. Even so, the total fractional amplitude 
can be as low as ~ l%-2% for a range of viewing angles 
only if the spots are within a few degrees of the stellar 
spin axis. For example, for our reference star the frac- 
tional modulation produced by two spots at the same 
inclination but separated by 160° in longitude is larger 
than the modulation produced by two antipodal spots if 
i > 70° and i sl = i s2 > 70° but is smaller if i < 60° 
and i s i = is2 ^ 60°. The total fractional amplitude can 
appear as low as ~ l%-2% over half the sky (i > 60°) 
only if the two spots are within ~ 10° of the spin axis. 

The intensity distribution produced by a Thomson 
scattering atmosphere is desc ribed by a linear c ombi- 
nation of Hopf functions (see iChadrasekhari 119601 Sec- 
tion 68) and peaks in the direction normal to the stellar 
surface. The radiation flux from such an atmosphere is 
given by the product of this intensity distribution and the 
projected area of the surface, which is proportional to the 
cosine of the angle to the normal. Consequently, the radi- 
ation flux from a Thomson scattering surface is strongly 
peaked in the direction normal to the surface. For this 
reason, emission from a Thomson scattering region typ- 
ically produces a higher modulation fraction than would 
isotropic emission from the same region. Hence if Thom- 
son scattering is important, the emitting regions must 
be closer to the star's spin axis than if the emission is 
isotropic. 

F an-shaped intensity dist r ibutio ns like those proposed 
by iPoutanen fc Gierlihskil (|2003f ) to describe Comp- 



tonized emission from near the stellar surface can par- 
tially compensate for the decrease of the projected area 
as the angle from the normal to the surface increases, pro- 
ducing a radiation flux distribution that is more isotropic 
than that produced by an isotropic intensity distribu- 
tion. For this reason, fan-shaped emission from a single 
spot can produce modulation at the first harmonic of the 
spin frequency that is smaller than would be produced by 
isotropic emission from the same spot. A spot produc- 
ing this emission pattern must still be close to the spin 
axis in order for the amplitude of the first harmonic to 
be < 2%. Fan-shaped emission from a single spot tends 
to increase the amplitude of the second harmonic rela- 
tive to the first. Fan-like emission from two antipodal 
spots produces waveforms similar to those produced by 
isotropic emission, but the modulation fraction can be 
up to a factor of 2 larger. 

In SAX J1808.4— 3658, the second harmonic is usually 
weak but is occasionally as strong as or even st ronger 
than the first harmonic (see lHartman et all [2008V The 
strengthening of the second harmonic could be caused 
by development of a two-fold asymmetry of one spot, 
by more fan-beam emission from a single spot, or by 
increased visibility a second, antipodal spot. 

Our computations show that modulation fractions are 
below typical current detection limits if the emitting re- 
gions are within 1° of the spin axis and remain there. 
As discussed in Section HI rapid fluctuations in the 
amplitude and phase of the oscillations would not be 
directly detectable in time-averaged observations, but 
would raise the background noise level and further re- 
duce the detectability of the oscillations. These effects, 
possibly in combination with others, such as reduction 
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Fig. 2. — Contour plots of the total fractional rms amplitude produced by isotropic emission from a single stable spot (left-hand panel) 
and from two stable antipodal spots (right-hand panel), as a function of the spot radius and the inclination of the primary spot relative to 
the spin axis, for a 1.4Mq star with a radius of 5M spinning at 400 Hz observed at an inclination of 45°. The heavier lines highlight the 
contours where the total flux modulation is 2%. This figure demonstrates that even very large spots must be centered close to the stellar 
spin axis in order to explain the low oscillation amplitudes observed in the AMXPs, unless almost the entire stellar surface is uniformly 
emitting. 



of the modulation fraction by scatter ing in circumstel- 
lar gas (|Lamb et al.lfl985t IMillenl2000h . may explain the 
nondetection of accretion-powered oscillations in accret- 
ing neutron stars in which nuclear-powered oscillations 
have been detected. 

3.2. Dependence on Spot Size 

Larger emitting spots tend to produce lower oscilla- 
tion amplitudes than smaller spots, but even very large 
spots must be centered close to the stellar spin axis in 
order to explain the low oscillation amplitudes observed 
in the AMXPs, unless almost the entire stellar surface 
is uniformly emitting. This is illustrated in Figure [21 
which shows the total fractional amplitude produced by 
isotropic emission from a single spot and from two an- 
tipodal spots, as a function of the radius of the spots and 
the inclination of the primary spot relative to the spin 
axis. All the curves in this figure are for our reference 
star, viewed at an inclination of 45°. 

The left-hand panel of FigureOshows that even a single 
spot that covers half of the star (a spot with a radius of 
90°) can produce a total fractional modulation as small 
as 2% only if its inclination is <5°. The modulation 
fraction depends only weakly on the size of the spot, 
especially for small-to- moderate spot inclinations. For 
example, the amplitude of the oscillation produced by 
a single spot inclined 45° from the spin axis decreases 
by only ~ 10% as the spot radius increases from 5° to 
45°. Even for a spot inclined 90° from the spin axis, 
the fractional flux modulation depends only weakly on 
the size of the spot for spot radii < 50°, decreasing from 
70% to 35% as the spot radius increases from 30° to 
90°. The reason for this insensitivity is that the angular 
width of the stellar radiation pattern produced even by 
spots with radii ~ 50° is governed mostly by the beaming 
pattern of the radiation from the stellar surface and the 
magnitude of the gravitational light deflection, not by 
the radius of the spot. A single spot at an inclination 
of 70° can produce a fractional modulation < 2% only if 
the entire surface of the star is uniformly emitting except 



for a dark area of radius 20° on the opposite side of the 
star from the center of the spot (i.e., if the spot radius is 
> 160°). 

The right-hand panel of Figure [5] shows that two uni- 
formly emitting antipodal spots of a given radius produce 
a smaller fractional modulation than a single spot of the 
same radius, because the antipodal spots cover twice as 
much of the stellar surface. Even so, the antipodal spots 
must be very large in order to produce fractional mod- 
ulations as small as those observed. For example, two 
antipodal spots on our reference star can produce a frac- 
tional modulation as low as ~ 2% for most observing di- 
rections only if they have radii > 75° , which means that 
all of the stellar surface is uniformly emitting except for a 
band around the star with a total width < 30° . The frac- 
tional modulation for most observing directions doubles 
to ^4% if the spot radius decreases from 75° to 65°. 

These results show that the low fractional amplitudes 
often observed in the AMXPs cannot be explained by 
large emitting regions, unless the emitting regions are 
uniform and cover almost the entire stellar surface, which 
is not expected. Another difficulty with attributing the 
low fractional amplitudes often observed to large emit- 
ting areas is that a substantial number of AMXPs that 
are observed to have fractional amplitudes ~ l%-2% at 
some times are observed to have much larger fractional 
amplitudes ~ 15%-20% at other times. Hence attribut- 
ing their small amplitudes at some times to large emit- 
ting areas would require their emitting areas to shrink 
by a large factor at other times, which is not expected. 

3.3. Dependence on Stellar Compactness 

The fractional modulation produced by a given emis- 
sion patt ern is generally smaller for more a compa ct star 
(see, e.g. IPechenick etal]|1983 IStrohmaverill992D . The 
reason is that gravitational focusing allows the observer 
to see a larger fraction of the surface of such a star, aver- 
aging the contributions of its bright and dark areas and 
reducing the flux variation seen as the star turns. Strong 
gravitational focusing can increase the fractional modu- 
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Fig. 3. — Contour plots of the total fractional rms amplitude produced by isotropic emission from a single stable spot (left-hand panel) 
and from two stable antipodal spots (right-hand panel), as a function of the stellar radius and the inclination of the primary spot relative 
to the spin axis, for a 1.4 Mq star with a radius of 5M spinning at 400 Hz observed at an inclination of 45°. The heavier lines highlight the 
contours where the fractional flux modulation is 2%. This figure demonstrates that compactness alone cannot explain the low-amplitude 
oscillations of the AMXPs. 



lation seen by some observers if the star is very compact 
(R<3.5M). 

Even if all AMXPs are relatively compact (R ~ AM), 
which is not expected, they would produce fractional 
modulations larger than the ~ l%-2% fractional mod- 
ulations often observed (see Section [J) unless their emit- 
ting regions are close to their stellar spin axes. Nor can 
the failure so far to detect accretion-powered oscillations 
in some nuclear-powered AMXPs be explained unless 
the accretion-powered emission comes from areas very 
close to the stellar spin axis. These points are illustrated 
by Figure [3l which shows the fractional amplitudes pro- 
duced by isotropic emission from a single spot (left-hand 
panel) and from two antipodal spots (right-hand panel) 
as a function of R/M and the inclination of the primary 
spot relative to the spin axis, for a I.AMq star spinning 
at 400 Hz observed at an inclination of 45°. 

The left-hand panel of Figure [3] shows that even if all 
the AMXPs had radii as small as 4M, their fractional 
modulations could be as small as ~ 2% for emission from 
a single spot only if the spot is within 3° of the spin 
axis. If instead the spot is at a moderate-to-high inclina- 
tion, their fractional modulations would be much larger. 
For example, if the spot is 45° from the spin axis, the 
fractional modulation would be ~ 25%. 

The right-hand panel of Figure [3] shows that although 
the fractional modulation produced by two antipodal 
spots is generally smaller than the modulation produced 
by a single spot, antipodal spots would have to be within 
10° of the spin axis to produce a fractional modulation 
as small as ~ 2%, even for a star as compact as AM. If 
instead the spots are 45°, from the spin axis, the frac- 
tional modulation would be ~ 7%, much larger than the 
lowest modulations observed in most AMXPs. 

As another example (not plotted here) , consider a rel- 
atively massive 2.2 M Q star with one or two spots that 
are not close to the spin axis, observed at an inclination 
of 45°. If there is a single, isotropically emitting spot 
of radius of 25° at an inclination of 45°, the fractional 
modulation would be ~ 30% if the star has a radius of 
4M and 22% if the star has a radius of 3.2Af. If instead 



there are two antipodal, isotropically emitting spots at 
inclinations of 45° and 135°, the fractional modulation 
would be ~ 11% if the star has a radius of AM and 9% 
if the star has a radius of 3.2M. These fractional mod- 
ulations are much larger than the lowest modulations 
observed in most AMXPs. The fractional modulation is 
typically even larger for stars that are still more com- 
pact, because of the strong gravitational focusing effect 
discussed above. 

These results show that even if the AMXPs are com- 
pact neutron stars, their fractional modulation will be 
much larger than the low fractional modulations they 
often produce, unless their emitting areas are within a 
few degrees of the spin axis. Hence the low fractional 
modulations often observed in the AMXPs cannot be at- 
tributed to high stellar compactness. 

A further difficulty in attributing the generally low 
fractional modulations of the AMXPs to high stellar 
compactness is that several of the AMXPs that exhibit 
fractional modulations ~ l%-2% at some times exhibit 
fractional modulations ~ 15%-25% only a few hours or 
days later. The compactness of the neutron star cannot 
change significantly on such short timescales, and hence 
some mechanism other than high compactness must be 
responsible for the low fractional modulations seen at 
many times in most AMXPs. 

Our results also show that high stellar compactness 
cannot be the main factor responsible for the nondc- 
tection of accretion-powered oscillations in many accret- 
ing neutron stars in low-mass X-ray binary systems, in- 
cluding many stars in which nuclear-powered oscillations 
have been detected. The upper limits on the fractional 
amplitude of any accretion-powered oscillation produced 
by these stars are < 0.5% (see Section[T]), much less than 
could be explained by the effects of high stellar compact- 
ness. 

3.4. Dependence on Stellar Spin Rate 

Most of the results discussed in this section are for stars 
spinning at 400 Hz. Other things being equal, stars spin- 
ning more rapidly will produce oscillations with larger 
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Fig. 4. — Fractional rms amplitudes of the full bolometric pulse profile (left) and the second harmonic (first overtone) component of 
the pulse profile (right) produced by isotropic emission from a single stable spot and from two stable antipodal spots, as functions of the 
inclination of the primary spot relative to the spin axis, for spot radii of 25° , a 1.4 Mq star with a radius of 5M spinning at 400 Hz, and the 
indicated observer inclinations i. This figure shows that if the inclination of the emitting region is small, modest changes in its inclination 
can produce large changes in the observed amplitude of the oscillation. 



fractional amplitudes, because their higher surface veloc- 
ities will produce larger Doppler boosts and greater aber- 
rati on, making their rad i ation patterns mo re asymmet- 
ric ([Miller fc LamrJll998t iBraje et alJl2000f ). Conversely, 
stars spinning more slowly tend to produce oscillations 
with smaller fractional amplitudes. However, the depen- 
dence of the fractional amplitude on the stellar spin rate 
is weak. Consequently, the basic conclusions reached in 
this section are valid for the full range of AMXP spin 
rates observed. 

4. PULSE AMPLITUDE AND PHASE VARIATIONS 

As discussed in Section [TJ the fractional amplitudes of 
the accretion-powered oscillations of the AMXPs vary by 
factors ranging from ~ 2 to ~ 10 on timescales of hours to 
days. As was also discussed in Section[TJ the phases of the 
Fourier components of these oscillations vary by ~0.T- 
0.4 cycles on similar timescales. If interpreted as caused 
by changes in the stellar spin rate, the observed phase 
variations would imply frequency variations more than a 
factor of 10 larger than expected for the largest accre- 
tion torques and the smallest inertial moments thought 
possible for these stars. In several AMXPs, the phase 
variations are correlated with the amplitude variations, 
for some amplitude ranges. A successful model of the 
accretion-powered oscillations of the AMXPs should pro- 
vide a consistent explanation of these properties of the 
oscillations. 

As noted in Section [1] changes in the latitude and 
longitude of the emitting area are expected on a wide 
range of timescales, from the ~ 0.1 ms dynamical time at 
the stellar surface to the ~ 10 d timescales of the varia- 
tions observed in the mass accretion rates of the AMXPs. 
Changes in the latitude (inclination) of the emitting area 
primarily affect the amplitudes of the harmonic compo- 
nents of the pulse profile but also affect their phases. 
Modest changes in the latitude of the emitting area can 
produce large changes in the oscillation amplitude if the 
emitting region is close to the stellar spin axis. Changes 
in the longitude (stellar azimuth) of the emitting area 
shift the phases of the harmonic components of the pulse 
but do not affect their amplitudes. 



In this section we show that if the emitting area is 
close to the spin axis, the fractional amplitude variations 
seen in all AMXPs and the larger fractional amplitudes 
~ 15%-20% occasionally seen in several of them can be 
explained by modest variations of the inclinations of their 
emitting areas. We show further that if the emitting area 
is close to the spin axis and moves in the azimuthal direc- 
tion by even a small distance, the phases of the Fourier 
components of the pulse will shift by a large amount. 
These shifts can easily be as large as ~ 0.1-0.4 cycles. In- 
deed, if the emitting area loops the spin axis, the phases 
of the Fourier components will shift by more than one cy- 
cle. We show that changes in the latitude and longitude 
of the emitting area by modest amounts can produce sub- 
stantial correlated changes in the amplitudes and phases 
of the Fourier components of the oscillation. 

4.1. Pulse Amplitude Variations 

Figure [4] shows the total fractional rms amplitude and 
the fractional rms amplitude of the second harmonic 
(first overtone) of the spin frequency for pulses produced 
by isotropic emission from a single stable spot and from 
two stable antipodal spots, as functions of the inclina- 
tion of the primary spot relative to the spin axis. These 
curves are for spots with angular radii of 25° on our ref- 
erence star. 

The results shown in Figure 3] support our conclusion 
in Section [3] that the oscillation amplitudes seen by most 
observers can be as small as ~ 1%— 2% only if the emitting 
regions are within a few degrees of the spin axis. If in- 
stead the emitting regions were > 45° from the spin axis, 
most observers would see modulation fractions > 15%. 

Whether the observer sees emission from a single area 
or from two antipodal areas, the total fractional am- 
plitude of the oscillation will increase if the inclination 
of an emitting region initially at a low inclination in- 
creases. The reason is that the star's radiation pattern 
becomes steadily more asymmetric about the spin axis 
as the emitting region moves farther from the spin axis, 
partly because the emitting region becomes more asym- 
metric and partly because the surface velocity becomes 
higher. 
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If the emitting area is initially close to the spin axis, 
a modest increase or decrease in its latitude can change 
the amplitude of the observed flux oscillation by a sub- 
stantial factor. For emission from a single spot, the am- 
plitude increase is steepest for observers at high inclina- 
tions, whereas for emission from two antipodal spots, the 
amplitude increase is steepest for observers at moderate 
inclinations. 

These effects are illustrated by the left-hand panel of 
Figure |H It shows that for a single spot initially at i a = 
10°, an increase of the spot inclination by as little as ~ 5° 
can increase the fractional amplitude of the oscillations 
seen by an observer at 45° by a factor of 2. For observers 
at higher inclinations, the oscillation amplitude increases 
even more. 

A single spot at a high inclination produces oscillations 
with a high fractional amplitude for observers at incli- 
nations > 60° and the amplitude increase only slightly 
as the inclination increases further. As an example, for 
the parameter values assumed in Figure [4] the total frac- 
tional amplitude seen by an observer at i = 45° increases 
steadily with increasing spot inclination, reaching ~ 70% 
when i s = 90°. In contrast, the fractional amplitude seen 
by an observer at i w 90° (i.e., close to the spin equator) 
is - 85% for i s = 70° and - 90% for i s = 90°. 

If instead the observer sees radiation from two antipo- 
dal spots, the total fractional amplitude increases with 
spot inclination for spot inclinations <50°, though not 
as steeply as for a single spot. The fractional ampli- 
tude eventually decreases with increasing spot inclina- 
tion, unless the observer is located in the spin equator. 
A change in the inclination of two antipodal spots by 
10° can change the total fractional amplitude seen by an 
observer at 45° by a factor of 2. 

The right-hand panel of Figure |4] shows that the frac- 
tional amplitude of the second harmonic is almost the 
same for radiation from one spot or from two antipo- 
dal spots, independent of the inclinations of the observer 
and the spot(s). Comparison of the right-hand panel of 
Figure H] with the left-hand panel shows that the second 
harmonic dominates the pulse profile for spot inclinations 
<80°. 

The accretion-powered oscillations observed in 
XTE J1814-338 and XTE J1807-294 occasionally 
have fractional amplitudes as large as 11% and 19%, 
respectively, althou gh they are us ually much s maller 
(iChung et a~l [2001 IPatnmol [200l iZhang et all [20061 : 
IChou et al.1 I2008fh Figure 2] shows that amplitudes 
this large are possible for isotropic emission from a 
single spot even for relatively small spot inclinations, 
provided the observer's line of sight is not too close 
to the star's spin axis. For example, amplitudes this 
large are possible for spot inclinations ~ 5°-15° if the 
observer's inclination is >50°. A decrease in the spot 
inclination to ~ l°-3° would then reduce the fractional 
amplitude to ~ l%-3%, the lowest values observed in 
these AMXPs. 

If the emitting areas of the AMXPs are typically lo- 
cated at high stellar latitudes, as in the model discussed 
here, their pulse amplitudes are expected to depend on 
the accretion rate and structure of the inner disk. These 
may vary on timescales as short or shorter than the 
~0.1 ms dynamical time at the stellar surface up to 
timescales as long as the ^ 10 d observed variations of 



the mass accretion rate. Consequently, the harmonic am- 
plitudes of AMXP pulses are expected to vary on these 
timescales. 

This model can be tested on timescales longer than 
the ~ 10 2 -10 3 s integration times required to construct 
pulse profiles by searching for correlations between the 
amplitudes of the harmonic components of the pulse 
profile or the variance of the pulse amplitude and the 
pulsar's X-ray flux or spectrum. Indeed, the rms frac- 
tional amplitudes of the pulses of XTE J1814— 338 ap- 
pear to be positivel y correlated with its mean X-ray flux 
(jChung et al.1 12008). Methods for testing the model on 
timescales shorter than the time required to construct a 
pulse profile are discussed in Section l4~5l 

4.2. Pulse Phase Variations 

A change in the longitude (stellar azimuth) of the emit- 
ting region alters the arrival time of the pulse. This shifts 
the measured phases of all the Fourier components by the 
same amount, relative to their phases if the pulse were 
produced by an emitting area fixed on the surface of a 
star rotating at a constant rate. Obviously, the phase 
of the pulse peak is also shifted by this amount. If the 
emitting area is close to the spin axis, a displacement 
in the azimuthal direction by even a small distance can 
produce large phase shifts. 

A change in the latitude (inclination) of the emitting 
area not only alters the amplitudes of the Fourier com- 
ponents of the pulse but also their arrival times (phases). 
These alterations change the shape and arrival time of 
the pulse in a complex way. They occur because a change 
in the latitude of the emitting area changes the velocity 
of the emitting surface, altering the aberration, Doppler 
shift, and propagation time of each photon. The result- 
ing changes in the pulse shape and arrival time depend 
on the beaming pattern of the radiation, the properties 
of the star, and the inclination of the observer. The 
phase shifts caused by changes in the latitude of the emit- 
ting area are expected to be modest, because the surface 
velocities of neutron star models with spin frequencies 
< 600 Hz are typically a small fraction of the speed of 
light. 

Figure O shows how the phases of the peak of the pulse 
and its first and second harmonic components vary as 
the emitting spot moves along a circular path on the 
stellar surface. Motion along this path changes both the 
latitude and the longitude of the spot. Motion in either 
direction along all or part of such a path may occur if the 
main impact area of the accreting matter moves around 
the magnetic pole as a result of changes in the structure 
of the inner disk, changes that could be caused by varia- 
tions of the accretion rate onto the star (see Section l2~Tj) . 

The left-hand panel of Figure \5\ introduces the an- 
gles used to describe the motion of the spot around 
the path. The right-hand panel shows the changes 
in the arrival times and phases of the peak of the 
pulse and its Fourier components that occur as the 
spot moves around the path shown in the left-hand 
panel. In this example, the phases of the peak and 
the first and second harmonic components vary by ~ 
0.35 cycles, which is comparab le to the phase y aria- 
tions seen in XTE J1807-294 (iMarkwardtl l2004f) and 
SAX J1808. 4-3658 (jMorgan et alj|2003t Hartma n et all 
2008) on timescales of hours to days. If the path of 
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Fig. 5. — Left: angles relevant to the path of the emitting spot considered here. The observer's inclination i is measured from the stellar 
spin axis Q. The path on the stellar surface along which the emitting spot moves is described in a coordinate system anchored in the star. 
The center of the circular path considered here is at colatitude 9 P and azimuth <f> p = in this coordinate system. The angular radius of the 
path is rj a . The azimuth £ s of the spot around the center of the path is measured from its lowest latitude. Right: phase and time residuals 
of the peak, first harmonic (fundamental), and second harmonic (first overtone) of the observed waveform as functions of the azimuth £ s 
of a circular spot wit h a r adius of 25° moving along a circular path like that shown in the left-hand panel. The phase and time residuals 
are defined in Section I2.3I In this example, the emission from the spot is assumed to be isotropic, p = 12° , r] a = 10° , i = 45° , the star 
has a mass of 1.4 Mq and a radius of 5M, and the star's spin rate is 400 Hz measured at infinity. 



the emitting region approaches or loops the spin axis, 
even a very small displacement can shift the phases by 
more than one cycle. The phase shifts shown in Figure [5] 
are typical for isotropic emission. The phase shifts for 
Comptoni zed-emission-beaming p attern s like those pro- 
posed by IPoutanen k. G icrliriski (2003)) are similar, but 
differ in detail. 

The phase and time residuals shown in Figure[5]are not 
zero when the spot is in the plane defined by the spin axis 
and the observer (£ s = 0) for two reasons. First, when 
£ s = the spot is not directly beneath the observer in 
this example, so a photon emitted at this moment will 
reach the observer later than would a photon emitted 
from the point directly beneath the observer, which is the 
emission point that defines the zero of the arrival time 
at the observer (see Section |2~3|) . Second, the spinning 
motion of the stellar surface produces aberrations, time 
delays, and Doppler energy shifts that distort the pulse, 
causing its harmonic components to reach the observer 
at times different from their arrival times from a nonro- 
tating star. In this particular example these effects cause 
the peak of the pulse and its harmonic components to ar- 
rive earlier than they would if the star were not rotating. 
The phase advance produced by these effects more than 
offsets the backward phase shift produced by the extra 
travel time. 

As the emitting spot moves around the path shown 
in Figure El its inclination increases from 2° to 22°. 
This causes the amplitude of the oscillation to in- 
crease from ~ 1% when the spot is closest to the rota- 
tion pole to ~ 15% when it is farthest from the pole. 
This variation is comparable to the amplitude varia- 
tions observed in X T E J1814-338 and XTE J1807-294 
dChung et all 12001 iPatrunol I2TM iZhang et all [200l 
IChou et al.H2008f ). 

As described above, the change in the inclination of the 
spot as it moves around the path shown in Figure [5] also 
contributes to the phase shifts. In general, an increase 
in the inclination of an emitting spot by 20° can by itself 



shift the phases of the Fourier components of the pulse by 
amounts ranging from ~ 0.05 to ~ 0.20 cycles, depending 
on the beaming pattern, stellar properties, and viewing 
direction, and can shift the phases of different harmonic 
components by different amounts. These effects are in- 
cluded in the phase shifts shown in Figure [5] 

These results show that if the emitting areas of the 
AMXPs are usually near their spin axes, as in the model 
discussed here, movement of the emitting area by mod- 
est distances can cause the arrival times of pulse Fourier 
components to vary by the large amounts observed, even 
if the AMXPs spin at a nearly constant rate. 

Emitting areas occur where accreting plasma impacts 
the stellar surface. As discussed in Section [27X1 the po- 
sition of these areas depends on how accreting plasma 
enters the magnetosphere and hence on the accretion 
rate and the structure of the inner disk. Consequently, 
pulse phase variations produced by changes in the lo- 
cation of the emitting areas should be correlated with 
changes in the luminosity and spectrum of the pulsar. 
The properties of the accretion flow are expected to vary 
on timescales ranging from as short as ~ 0.1 ms to as long 
as ~ 10 d, and the arrival times of pulse Fourier compo- 
nents should therefore vary on these timescales also. 

This model of pulse phase variations can be tested 
on timescales longer than the ~ 10 2 -10 3 s times re- 
quired to construct a pulse profile by searching for and 
studying correlations between the measured phases of 
the harmonic components of pulses or the variance of 
these phases and the pulsar's X-ray flux or spectrum. 
The phase residuals of the first and second harmonic 
components of the pulses of XTE J1814 — 338 appear 
to be anticorrelated with its X-ray flux (|Chung et al.1 
120081 iPapitto etTaTl 120071 ) , making it a good candidate 
for this kind of study. iRiggio et alJ (2008) have re- 
ported a strong anticorrelation between the X-ray flux 
of XTE J1807-294 and the phase residual of the sec- 
ond harmonic component of its pulse profile. Methods 
for testing this model of phase variations on timescales 
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Fig. 6. — Distribution of pulse arrival times as a function of the 
pulse amplitude expected for pulses produced by an emitting area 
that is close to the star's spin axis and wanders. In this example 
the emitting area is assumed to wander in stellar longitude by the 
same distance di at all latitudes. The dotted and solid curves show 
the expected scatter of pulse arrival times (right vertical axis) and 
arrival phases (left vertical axis) as a function of the pulse ampli- 
tude, assuming / = di/2-irR = ±0.001. The dashed curve shows 
the relation between the arrival time and amplitude of pulses pro- 
duced by a spot that remains at the same longitude as its latitude 
changes. Here the arrival time of a pulse is taken to be the arrival 
time of i ts peak. The phase and time residuals arc defined in Sec- 
tion 12.31 These results are for a spot with a radius of 25° on a 
1.4 Mq star with a radius of 5M spinning at 400 Hz, observed at 
an inclination of 45° . The amplitude range plotted on the horizon- 
tal axis corresponds to spot inclinations from 2° to 45° (compare 
Figure [2). This figure illustrates the large scatter expected when 
the fractional amplitude is small. 

shorter than the time required to construct a pulse profile 
are discussed in Section l4~5l 

4.3. Correlated Pulse Amplitude and Phase Variations 

If the emitting areas of the AMXPs are near their spin 
axes and move around with time, the amplitudes and 
phases of their pulses should show two types of correlated 
behavior. 

A strong expectation in this model is that the scatter 
of the pulse arrival times (i.e., the scatter of the time or 
phase residuals) should decrease steeply with increasing 
pulse amplitude. The reason for this is that, as shown 
previously, a given variation in the azimuthal position 
of the emitting area produces a much larger phase shift 
when the emitting area is very close to the spin axis 
than when it is farther away; in contrast, the oscillation 
amplitude is much smaller when the emitting area is very 
close to the spin axis than when it is farther away. 

This effect is illustrated in Figure which shows the 
expected distribution of pulses in the arrival time versus 
amplitude plane if they are produced by emission from 
an area close to the star's spin axis that wanders over 
time. In the example shown, the phase residuals vary by 
~ 0.3 cycles when the fractional amplitude is ~ 0.02 but 
by only ~ 0.03 cycles when the fractional amplitude is 
~0.1. The pulse phase resi duals seen in s everal AMXPs 
(see, e .g., XTE J1807 -294: |Patruncj[2008; SAX J1748.9- 
2021 : IPatrund[200l and IGR J00291+5934: IPatrunol 
2008) have distributions similar to that shown in Fig- 
ure [51 These distributions are consistent with the model 
discussed here if the emitting areas of these AMXPs 
wander in the azimuthal direction by distances equal to 



~ 10~ 2 -10~ 3 of the stellar circumference. 

A second expectation in the moving spot model dis- 
cussed here is that the pulse arrival time or phase resid- 
uals are likely to form a track in the phase-residual ver- 
sus pulse-amplitude plane, especially if the change in the 
pulse amplitude is large. Such a track will be formed if 
the emitting areas where the accreting matter impacts 
the stellar surface tend to move along a favored path in 
stellar latitude and longitude as the accretion rate and 
the structure of the inner disk change. This is expected 
because models of the flow of accreting matter from the 
inner disk to the stellar surface predict that matter will 
be guided along particular but different magnetic flux 
tubes as the accretion rate and the structure of the inner 
disk vary (see Section 12. 1[) . 

Our computations show that motion of the emitting 
area along a path that changes its latitude but not its 
longitude can shift the phases of the first and second har- 
monics by at least 0.15 cycles and by different amounts. 
Motion of the emitting area along a path that instead 
changes its longitude but not its latitude can shift the 
phases of the first and second harmonics and the pulse 
peak by much larger amounts, if the emitting area is near 
the spin axis, but in this case the phase shifts will be the 
same for all Fourier components. Only if the shift in 
pulse phase produced by the change in the longitude of 
the emitting area exactly compensates for the shift pro- 
duced by the change in the latitude of the emitting area 
will there be no trend in the pulse phase residuals with 
increasing pulse amplitude. 

No track may be discernible in the phase-residual ver- 
sus pulse-amplitude plane when the pulse amplitude is 
small because, as has just been discussed, the pulse phase 
is very sensitive to the position of the emitting area when 
the area is near the spin axis, which is where it is ex- 
pected to be when the pulse amplitude is small. If, how- 
ever, the emitting area moves well away from the spin 
axis, which in the moving spot model produces a sub- 
stantial increase in the pulse amplitude, a track may be- 
come apparent. Tracks are observed in the dat a on sev- 
eral AMXPs when the pulse amplitude is > 3% (jPatrunol 
2008). The phase of the first harmonic component of 
the pulse is correlated with the amplitude of the pulse 
in IGR J00291+5934 and XTE J1751-305 but anticor- 
related in XTE J1807-294 and XTE J1814-338. 

The path on the stellar surface along which the emit- 
ting area moves as the accretion rate changes depends 
on the details of the accretion flow from the inner disk 
to the stellar surface that cannot yet be determined from 
first principles. However, the path along which the emit- 
ting area moves in the moving spot model can be inferred 
from simultaneous measurements of the phases and am- 
plitudes of pulses. 

Similar shifts in the phases of all Fourier components 
are consistent with movement of an emitting area of 
roughly fixed size, shape, and radiation-beaming pat- 
tern along a path that causes a substantial change in 
its longitude as its latitude decre ases. This behav- 
ior is observed in XTE J1814-33 8 (jPapitto et all 120071 : 
iChung et all 120081 : IPatrunol l200l ). If instead the phase 
residuals of different Fourier components of the pulse 
evolve very differently with increasing pulse amplitude, 
this is an indication that the emitting area is moving 
along a path that produces very little change in the area's 
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longitude and/or a substantial change in its shape or the 
radiation-beamin g pattern . This behavior is observed in 
XTE J1807-294 |Patrundl20"08h . 

4.4. Accretion- and Nuclear-Powered Oscillations 

The close agreement of the pulse profiles and phases of 
the accretion- and nuclear-powered ( X-ray burst) oscilla- 
tions observed in SAX J1808.4 -3658 (iChakrabartv et all 
120031 ) and XTE J1814-338 (jStrohmaver et alJ 12003] ) 
strongly suggests that in these stars, both types of os- 
cillation are produced by X-ray emission from nearly the 
same area on the stellar surface. If this is so, it implies 
that in these AMXPs thermonuclear burning is concen- 
trated near the magnetic poles onto which accreting mat- 
ter is falling. It also implies that long-term variations in 
the phase residuals of the two types of oscillations should 
track one another in these pulsars and should also be cor- 
related with variations in the X-ray flux and spectrum, 
because both types of variations are produced by changes 
in the accretion flow through the inner disk. 

This interpretation can be tested by comparing the 
observed variations of the amplitudes and phases of 
the Fourier components of the accretion- and nuclear- 
powered oscillations with the variations predicted by this 
model. If this interpretation proves correct, the locations 
and movements of the emitting areas can be determined 
from the observed phase and amplitude variations. 

In Section 15.11 we emphasize that a mechanism that 
drives a star's magnetic poles toward its spin axis can 
greatly reduce the dipole component of the star's mag- 
netic field without reducing significantly its strength. 
The final configuration of the magnetic field produced 
by this mechanism could have a total strength ~ 10 11 - 
10 12 G, but a dipole component ~ 10 s G, consistent with 
the dipole moments inferred from the X-ray emission and 
spin evolution of AMXPs and millisecond radio pulsars 
(jLamb fe Boutloukosl [20081 ). 

Surface magnetic fields ~ 10 11 -10 12 G are strong 
enough to laterally confine accreting matter within the 
surface layers of the neutron star, creating an accu- 
mulation point for fuel for thermonuclear bursts and 
causing thermonuclear bur ning to be concentrated near 
the star's magnetic poles dWooslev fe Wallace! 1 1982t see 
also iMelatos fe Phinnevl 120011 and references therein) . 
They are also strong enough to produce strong-magnetic- 
field features in the keV X-ray spectra of AMXPs (see 
iMeszarosl 119921 ) . If this picture is correct, such fea- 
tures are more likely to be detected in AMXPs such as 
SAX J1808.4-3658 and XTE J1814-338 that produce 
nuclear-powered oscillations that are phase-synchronized 
(or nearly synchronized) with their accretion-powered os- 
cillations than in other AMXPs. 

4.5. Effects of Rapid Spot Movements 

Our results show that motion of the emitting area on 
the stellar surface generally changes both the amplitudes 
and the phases of the Fourier components of the pulse 
profile. As noted in Section |2~T1 the position of the emit- 
ting area is expected to reflect the accretion rate and 
structure of the inner disk, and is therefore expected to 
change on timescales at least as short as the ~0.1 ms 
dynamical time near the neutron star. Pulse phase and 
shape variations are therefore expected on this timescale. 



The observable effects of changes in the position of the 
emitting area on timescales longer than the ~ 10 2 -10 3 s 
integration times required to construct a pulse profile 
have been discussed in previous sections. Here we dis- 
cuss the expected effects of fluctuations in the position 
of the emitting area on timescales shorter than the time 
required to construct a pulse profile. 

Changes in the position of the emitting area on 
timescales shorter than the time required to construct 
a pulse profile w ill appear as unre solved amplitude and 
phase noise (see lLamb et al.l 119851 ) . This noise will re- 
duce the apparent fractional amplitudes of the oscilla- 
tions at the spin frequency and its harmonics, making 
them appear weaker than in the plots presented in this 
paper, which assume stable emitting spots fixed on the 
stellar surface. These pulse phase and amplitude fluc- 
tuations are likely to be greater when the emitting area 
is near the spin axis, because there a displacement by 
a given distance produces a larger change in the pulse 
phase and, for many geometries, in the pulse amplitude. 
This trend will tend to reduce further the apparent pulse 
amplitude when the emitting area is close to the spin 
axis, making it appear even smaller relative to the pulse 
amplitude when the emitting area is further from the 
spin axis than is shown in the plots presented earlier in 
this section. 

Pulse shape fluctuations caused by rapid changes in 
the position of the emitting region will appear as back- 
ground noise in excess of the Poisson counting noise. 
This noise can in principle be detected, especially be- 
cause its strength is expected to be anticorrelated with 
the pulse amplitude and to depend in a systematic way 
on the X-ray flux and spectrum of the pulsar. Detection 
of more excess noise when the apparent pulse amplitude 
is smaller would support the moving spot model of pulse 
phase and amplitude variations. 

4.6. Undetected and Intermittent Pulsations 

The results presented in Section [3] show that if the 
emitting area is very close to the spin axis and re- 
mains there, the amplitude of X-ray oscillations at the 
stellar spin frequency or its overtones may be so low 
that they are undetectable. In addition, rapid vari- 
ations in the shape and phase of pulses are expected 
to be stronger when the emitting area is very close to 
the spin axis. The noise produced by these fluctuations 
may — in combination with other effects, such as reduc- 
tion of the modulation fraction by scattering in circum- 
stellar gas — further reduce the detectability of accretion- 
powered osci llations in neutron stars with millisecond 
spin periods (jLamb et al.lll985l : iMillerl 12000ft . Taken to- 
gether, these effects may help explain the nondctection of 
accretion-powered oscillations in some accreting neutron 
stars in which nuclear-powered oscillations have been de- 
tected. 

Location of the X-ray emitting areas close to the spin 
axis also suggests a natural explanation for the sudden 
appearance and disappearance of accretion-powered os- 
cillations observed in the "intermittent" AMXPs (see 
lLamb et alJ I2008at [Boutloukos et all 120081 ). A change 
in the accretion flow within the magnetosphere can sud- 
denly channel gas to the stellar surface farther from the 
spin axis, causing the centroid of the emitting area to 
move away from the axis. This will increase the ampli- 
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tude of the oscillation, potentially making a previously 
undetectable oscillation detectable. Conversely a change 
in the accretion flow that suddenly channels accreting gas 
closer to the spin axis could make a detec t able o scilla- 
tion undetectable. As noted bv lLamb et al.l (|2008af ). this 
idea can be tested by studying short-term variations in 
the amplitudes and phases of the harmonic components 
of the pulses. 

5. DISCUSSION 

The results presented in previous sections show that a 
model of AMXPs in which the X-ray-emitting areas are 
close to the spin axis and move around on the stellar sur- 
face can explain many of their properties. Emitting areas 
close to the spin axis are to be expected if the magnetic 
poles of the AMXPs are close to their spin axes, caus- 
ing accreting gas to be channeled there. In this section, 
we first discuss mechanisms that may cause the mag- 
netic poles of AMXPs to be close to their spin axes. We 
then point out that this picture of magnetic field evolu- 
tion may also explain why AMXPs in which accretion- 
powered oscillations have been detected are transient X- 
ray sources. We also discuss several observational tests 
of this model. 

5.1. Movement of Magnetic Poles Toward the Spin Axis 

The neutron vortices in the fluid core of a spinning 
neutron star are expected to move radially inward if the 
star is spun up. This inward vortex motion is expected 
to drag the magnetic flux tubes in the fluid core of the 
star toward the star's spin axis (|Srinivasan et al.lll99fl ; 
I Ruder man! 1 1 99 If ) . 4 Some of the relevant physics is not 
well understood and many details remain to be explored, 
but this process is expected to squeeze the magnetic flux 
of accreting pulsars toward their spin axes as they are 
spun up to short periods (recycled). As a result, the mag- 
netic poles of recycled millisecond p ulsars are expected 
to be very close to the ir spin axes (jChen fc Rudermanl 
119931: iChen et al.lll998t) . The magnetic poles of recy- 
cled pulsars could also end up close to their spin axes if 
the accretion spin-up torque has a component that tends 
to align the star's magnetic field with its spin axis as 
it spins up. Motion of the magnetic poles toward the 
spin axis will be facilitated by diffusion of magnetic flux 
through t he crust and/or f racture and fragmentation of 
the crust (|Rudermanl 1 1 99 If ) . 

If the star's north and south magnetic poles are in op- 
posite rotation hemispheres when spin-up begins, the in- 
ward motion of vortices will drag them toward opposite 
spin poles. If instead both poles are in the same rotation 
hemisphere when spin-up begins, they will be dragged 
toward the same spin pole. In either case, the strength 
of the magnetic field at the stellar surface will be orders 
of magnitude larger than the strength of its dipole com- 
ponent, as we now explain. 

An initial magnetic field M\ that is approximately uni- 
form over a stellar hemisphere of radius R has a magnetic 
moment \i\ w M%R '. The magnetic flux threading the 

4 It has also been argued (Srinivasan et al. 1990) that a neutron 
star's magnetic field will be greatly reduced earlier in its evolution if 
it is spun down by a large factor, causing outward moving neutron 
vortices to force its magnetic field into the crust near the rotation 
equator, where it can be dissipated. 



fluid core is conserved as the magnetic field is squeezed. 
Suppose the inward moving vortices squeeze the flux into 
a tube of radius a « R/30. This will produce a squeezed 
magnetic field of strength M 2 » (R/a) 2 Mi » 10 3 M X . 

If the star's north and south magnetic poles are forced 
toward opposite spin poles, the dipole moment of the 
squeezed magnetic field will be aligned with the star's 
spin axis and will have a magnitude fi2 ~ M2R0? ~ fix, 
corresponding to a surface magnetic field with a dipole 
component M2 (dipole) w Mx- The full strength M2 of 
the surface magnetic field will therefore be w (R/a) 2 w 
10 3 times larger than the strength M2 (dipole) of its 
dipole component. 

If instead the star's north and south magnetic poles 
are both forced toward the same spin pole, the dipole 
moment of the squeezed magnetic field will be orthog- 
onal to the spin axis and hence the star will be an 
"orthogonal rotator" , even though both magnetic poles 
are near a spin pole rather than near the spin equa- 
tor. In this case the dipole moment will have a mag- 
nitude fi2 ~ M2L0? w [L/R)fix, where L is the final 
separation between the north and south magnetic poles, 
corresponding to a surface magnetic field with a dipole 
component M 2 (dipole) w (L/R)Mx. Thus if L w 0.1R, 
the dipole moment of the squeezed field will be about 
10 times smaller than the dipole moment of the original 
field and the full strength M2 of the surface magnetic 
field will be « (R/L)(R/a) 2 « 10 4 times larger than the 
strength M 2 (dipole) of its dipole component. 

Regardless of whether the star's north and south mag- 
netic poles are forced toward opposite spin poles, cre- 
ating an aligned rotator, or toward the same spin pole, 
creating an orthogonal rotator, both magnetic poles will 
end up very close to the spin axis. Accreting gas that is 
channeled along magnetic field lines will therefore impact 
the stellar surface close to the spin axis and accretion- 
powered X-ray emission will come primarily from regions 
near the spin axis. 

Note that if a pulsar's north and south magnetic poles 
have both been forced close to the same spin axis, the 
flow of accreting matter could impact the surface in a 
pattern with predominantly either one-fold or two-fold 
symmetry relative to the spin axis. In this accretion 
geometry, even a small change in the properties of the 
flow through the inner disk could change the dominant 
symmetry of the emitting area from one-fold to two-fold 
or vice versa, thereby changing the relative strengths of 
the first and second harmonics in the pulse profile. 

Another consequence of this picture of magnetic field 
evolution is that both accretion-powered (X-ray) and 
rotation-powered (radio) millisecond pulsars may have 
surface magnetic fields ~ 10 n -10 12 G. This is - 10 3 -10 4 
times stronger than the ~ 10 8 -10 9 G dipole components 
infer red from the magnetospheric r adii of AMXPs (see, 
e.g., iPsaltis fc Chakrabartvl Il999t lLamb fc Boutl oukos 
EOQi), the spin-down of SA X J1808.4-3658 in quies- 
cence (|Hartman et al.l |2008|) , and the spin-down rates 
of the rotation-powered m illisecond radio pulsars (sec 
lLamb fc Boutloukosll2008r ). Magnetic fields this strong 
are strong enough to create an accumulation point 
for fuel for thermonuclear bursts and produce strong- 
magnetic-field features in the keV X-ray spectra of 
AMXPs (see Section |4~4]) . 
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5.2. Why AMXPs are Transients 

The picture of the X-ray emission and magnetic field 
evolution of AMXPs that we have outlined here suggests 
a possible explanation for why all AMXPs found so far 
are transients. These systems have very low long-term 
average mass transfer rates, but binary modeling sug- 
gests that the mass transfer rates were h igher in the past 
(see, e.g.. iBildsten fc Chakra bartv 2001 for a discussion 
in the context of SAX J1808.4-3658). If stars such as 
these are initially spun up to high spin rates, so that their 
magnetic poles are forced very close to their spin axes, 
they would appear similar to the accreting neutron stars 
in low-mass X-ray binary systems in which accretion- 
powered oscillations have not been detected. When their 
accretion rates later decrease, and magnetic dipole and 
other braking torques cause them to spin down, their 
magnetic poles will be forced away from the rotation axis, 
and their accretion-powered oscillations will become de- 
tectable. 

If this explanation of the transient nature of the 
AMXPs is correct, they should be spinning down on 
long timescales, as seems to be the case (see, e.g., 
lHartman et al.|[2008h . and the amplitudes and phases the 
harmonic components of their waveforms should evolve 
in a correlated way as their magnetic poles move away 
from their spin axes (see Section HJ. 

5.3. Comparison with the Properties of MRPs 

Rotation-powered millisecond radio pulsars (MRPs) 
are thought to be the offspring of AMXPs and therefore 
should have magnetic field geometries similar to those 
of the AMXPs, except that the low rate of accretion at 
the end of the accretion phase and magnetic dipole brak- 
ing toward the end of the accretion phase and afterward 
will tend to spin them back down, forcing their magnetic 
poles away from their spin axes. 

Although the radio emission properties of most MRPs 
provide little clear evidence about the locatio ns of their 
magnetic poles fsee iManchester fc Hani 120041 ) . there are 
indications from their radio emission that the dipole mo- 
ments of a substantial number of the most rapidly spin- 
ning MRPs in the galactic disk are nearly aligned or 
nearly orthogonal to their spin a xes (|Chen et al.l 119981 : 
see also lChen fc Rudermadfl993l ). As explained above, 
either orientation is consistent with their magnetic poles 
being driven close to their spin axes by neutron vortex 
motion during spin-up as AMXPs. 

Analysis and modeling of the waveforms of the ther- 
mal X-ray emission from MRPs may provide better con- 
straints on their magnetic field geometries. Recent mod- 
eling of the high (30% to 50% rms) amplitude X-ray os- 
cillations observed in three nearby MRPs using Comp- 
tonized emission and two antipodal or nearly antipodal 
hot spots is consistent with their emitting regions being 
far from their spin axes (jBogdanov et al.ll2007ll2008f) . We 
note that all three pulsars have relatively low (~ 150 Hz- 
200 Hz) spin rates and may therefore have been spun 
down by a factor ~ 3 from their maximum spin frequen- 
cies after spin-up. If so, their magnetic poles may have 
been forced away from their spin axes by a similar fac- 
tor. This could allow co nsistency between the ~ 40°-60° 
inclinations inferred by iBogdanov et al.l for these three 
MRPs and the < 20° inclinations typically required in 



the model of AMXPs described here. We also note that 
only MRPs that produce high-amplitude X-ray oscilla- 
tions are currently detectable as oscillators, so the ones 
that are detected may have magnetic inclinations larger 
than is typical. More investigation is needed. 

5.4. Other Observational Tests 

In the nearly aligned moving spot model of AMXP X- 
ray emission, the properties of X-ray pulses (e.g., their 
amplitudes, harmonic content, and arrival times) should 
be functions of the pulsar's X-ray luminosity and spec- 
trum. The reason is that the properties of pulses are 
determined by the location of the emitting area, which 
depends on the rate and structure of the accretion flow 
through the inner disk. The X-ray luminosity and spec- 
trum of the pulsar will also depend on the accretion flow 
through the inner disk. Hence the properties of pulses 
should be correlated with the X-ray luminosity and spec- 
trum. 

Changes in the location of the emitting area on the 
stellar surface on timescales longer than the ~ 10 2 -10 3 s 
intervals required to construct stable pulse profiles will 
produce correlated changes in the amplitudes, harmonic 
content, and arrival times of these profiles. 

In the nearly aligned moving spot model, the arrival 
times of pulses with low amplitudes are expected to fluc- 
tuate much more than the arrival times of pulses with 
high amplitudes (see Figure [5]) . The reason is that pulse 
amplitudes are lower when the emitting area is closer to 
the spin axis, where small variations in the position of the 
emitting area can produce large variations in the pulse 
arrival time. 

A second expectation is that pulse phase residuals will 
form a track in the phase-residual versus pulse-amplitude 
plane. This will be the case if the area where accreting 
matter impacts the stellar surface and X-rays are emit- 
ted tends to move along a favored path on the stellar sur- 
face as the accretion flow through the inner disk changes. 
The location of the emitting area on its favored path 
will depend on the accretion flow through the inner disk 
and will in turn determine the position of pulses in the 
phase-residual versus pulse-amplitude plane. Changes in 
the accretion flow through the inner disk will therefore 
cause correlated changes in the phases and amplitudes 
of pulses, producing a track when these are plotted in 
the phase-residual versus pulse-amplitude plane. Such 
tracks are likely to be more evident if the range of the 
pulse amplitude variation is large. 

The pulses of IGR J00291+5934, XTE J1751-305, 
XTE J1807-294, and XTE J1814-338 seem to form 
track-like patterns in the phase-residual versus pulse- 
amplitude plane when the pulse amplitude is > 5% (see 
Section |4.3[) . An important expectation of the moving 
spot model is that the position of pulses along such a 
track should be correlated with the X-ray luminosity and 
spectrum of the pulsar. 

Another check of the moving spot model is possi- 
ble if the nearly identical pulse profiles and phases of 
the accretion- and nuclear-powered X-ray oscillations of 
XTE J1814-338 and SAX J1808.4-3658 are produced 
by emission from the same area on the stellar surface, as 
suggested in Section 14.41 If this is the case, the model 
predicts that longer-term (days-weeks) variations of the 
phase residuals of the accretion- and nuclear-powered os- 
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cillations should be correlated with one another and with 
longer-term variations of the pulsar's X-ray luminosity 
and spectrum, because all three variations are expected 
to be related to changes in the accretion flow through 
the inner disk. Evidence of this behavior would support 
the moving spot model. 

The location of the emitting area on the stellar surface 
is likely to change on timescales as short as the ~ 1 ms 
dynamical timescale in the inner disk. These rapid move- 
ments of the emitting area will cause rapid variations of 
the observed X-ray flux. Flux variations on timescales 
shorter than the ~ 10 2 -10 3 s intervals required to con- 
struct stable pulse profiles are observable as noise in ex- 
cess of the photon counting noise. In the moving spot 
model, the motion of the emitting area that produces 
this excess noise also affects the amplitudes, harmonic 
content, and arrival times of the X-ray pulses, and is 
related to the X-ray luminosity and spectrum. Hence 
the strength of the excess noise produced by motion of 
the emitting area should be correlated with these other 
properties of the AMXP. 

A strong expectation in the nearly aligned moving spot 
model is that the component of the excess noise produced 
by motion of the emitting area should be stronger when 
the pulse amplitude is smaller and weaker when the pulse 
amplitude is larger (see Section rO|) . 

One way to focus on the strength of the noise pro- 
duced by movement of the emitting area would be to fil- 
ter the X-ray flux data to remove variations on timescales 
shorter than, say, 0.05 s. This would remove the vari- 
ability associated with oscillations at the pulsar spin 
frequency and any high-frequency QPOs. The remain- 
ing total variability would then be correlated with the 
other properties of the AMXP listed above. Detection 
of greater variability when the observed pulse amplitude 
is smaller would support the nearly aligned moving spot 
model. 

If AMXPs do have surface magnetic fields as strong 
as ~ 10 11 -10 12 G, as suggested in Section 15.11 their 
keV spectra may show strong-magnetic-field features. 
Such features may include increased flux below the fun- 
damental cyclotron resonance energy, where the mag- 
netic fiel_d_sj : n3£resses_the opacity of surface layers (see, 
e.g., iPavlov et alj fl975): resonance scattering features 
produced by the funda mental cyclotron resonanc e and 
its overtones (see, e.g., IMeszaros fc Naeell fl985h : and 
changes in the spectrum at the higher energies where 
the cyclotron scattering becomes unimportant. 

Another implication of the existence of AMXPs with 
surface magnetic fields as strong as ~ 10 11 -10 12 G is that 
the rotation-powered millisecond pulsars that are their 
progeny should have surface magnetic fields this strong, 
even though their dipole fields are ~ 10 8 -10 9 G. Mag- 
netic fields this strong can be expected to affect particle 
acceleration and 7-ray emission by these pulsars. 

6. CONCLUSIONS 

In previous sections we have explored in some detail 
the nearly aligned moving spot model of AMXP X-ray 
emission. In this model the X-ray emitting regions are 
close to the stellar spin axis and move around on the 
stellar surface with time. Here we list our principal con- 
clusions. 

Pulse amplitudes and shapes. In Section 13. 1 1 we in- 



vestigated the amplitudes and shapes of the pulses pro- 
duced by emitting spots on the stellar surface as a func- 
tion of their inclination to the spin axis, for several X-ray 
beaming patterns and a range of stellar masses, compact- 
nesses, and spin rates. We found that emitting areas on 
or near the stellar surface can produce fractional ampli- 
tudes as low as the l%-2% values often observed only if 
they are located within a few degrees of the stellar spin 
axis. Regions near the spin axis also naturally produce 
nearly sinusoidal pulse profiles. 

We explored effect of spot size on pulse amplitude 
in Section 13.21 We found that although the pulse am- 
plitudes produced by large emitting areas tend to be 
smaller, this effect is weak. Unless almost the entire sur- 
face of the star is uniformly emitting, even large spots 
produce pulse amplitudes greater than those observed in 
the AMXPs, unless the spots are centered close to the 
spin axis. 

In Section 13 . 31 we studied the effect of stellar com- 
pactness on pulse amplitudes. Although the pulse am- 
plitudes produced by very compact neutron stars tend to 
be smaller than the amplitudes produced by less compact 
stars, we found that this effect is too weak to explain by 
itself pulse amplitudes as small as those observed in the 
AMXPs. Stellar compactness clearly cannot explain why 
the fractional amplitudes of several AMXPs are ~ 1%- 
2% at some times but ~ 15%-25% a few hours or days 
later, because the stellar compactness cannot change on 
such short timescales. 

These results show that emission from the stellar sur- 
face can explain the low amplitudes and nearly sinusoidal 
waveforms typically observed in AMXPs only if the emit- 
ting areas are located close to the stellar spin axis. 

Variability of pulse amplitudes, shapes, and arrival 
times. In Section 14.11 we investigated the amplitude 
changes that can be produced by motion of the emitting 
area on the stellar surface. We found that if the emit- 
ting area is close to the spin axis, even a small change 
in the latitude of the area can change the oscillation am- 
plitude by a substantial factor. For example, changes 
in the inclination of the emitting area by < 10° can ex- 
plain the amplitude variations seen in the AMXPs and 
the relatively large fractional amplitudes ~ 15%-20% oc- 
casionally seen in some of them. 

In Section 14. 2i we studied the changes in the arrival 
times (phases) of the harmonic components of the pulse 
caused when the emitting area moves around on the stel- 
lar surface. We found that changes in the latitude of the 
emitting area can shift the phases of the first and sec- 
ond harmonics by at least 0.15 cycles and by different 
amounts. We showed that if the emitting area is close 
to the spin axis and moves in the azimuthal direction by 
even a small distance, the phases of the first and second 
harmonics can easily shift by as much as ~ 0.1-0.4 cycles. 
If the emitting area loops the spin axis, the phases of the 
Fourier components will shift by more than one cycle. 

Motion of the emitting area on the stellar surface gen- 
erally produces both amplitude and phase variations. As 
discussed in Section [21 the position of the emitting area 
is expected to reflect the accretion rate and structure 
of the inner disk, and is therefore expected to change 
on timescales at least as short as the ~0.1 ms dynami- 
cal time at the stellar surface and as long as the ~ 10 d 
timescale of the variations observed in the mass accretion 
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rate. Changes in the position of the emitting area on 
timescales longer than the ~ 10 2 -10 3 s integration times 
required to construct a pulse profile will produce changes 
in the apparent pulse amplitude and phase. 

These results show that if the emitting area is close to 
the spin axis, modest changes in its location can explain 
the rapidly varying harmonic amplitudes and phases of 
the AMXPs. 

Correlated amplitude and phase variations. We showed 
in Section 14.31 that changes in the latitude and lon- 
gitude of the emitting area tend to produce corre- 
lated changes in the amplitudes and phases of the har- 
monic components of the pulse. A strong expecta- 
tion in the nearly aligned moving spot model is that 
the scatter in the pulse arrival times (i.e., the pulse 
time or phase residuals) should decrease steeply with 
increasing pulse amplitude. The residuals of several 
AMXPs, including XTE J1807-294, SAX J1748.9-2021, 
and IGR J00291+5934, behave in this way. The mag- 
nitudes of the phase residuals of these AMXPs are con- 
sistent with the nearly aligned moving spot model if the 
emitting areas wander in the azimuthal direction by dis- 
tances ~ 10 -2 -10 -3 times the stellar circumference. 

A second expectation in the nearly aligned moving spot 
model is that the arrival times of pulses will form a track 
in the phase-residual versus pulse-amplitude plane, es- 
pecially if the change in the pulse amplitude is large. 
Such a track will be formed if the emitting areas where 
the accreting matter impacts the stellar surface move re- 
peatedly along a particular path in stellar latitude and 
longitude as the accretion rate and the structure of the 
inner disk change. This is expected because models of the 
flow of accreting matter from the inner disk to the stellar 
surface predict that matter will be guided along partic- 
ular but different magnetic flux tubes as the accretion 
rate and the structure of the inner disk vary (see Sec- 
tion [23}. A s discussed in Section |4~31 tracks of this type 
are observed in plots of pulse arrival time versus pulse 
amplitude for XTE J1814-338 and XTE J1807-294. 

Accretion- and nuclear-powered oscillations. The suc- 
cess of the emitting spot model discussed here supports 
magnetic field evolution models in which the magnetic 
flux of the accreting neutron star becomes concentrated 
near its s pin a xis as it is spun up. As discussed in Sec- 
tions 14.41 and 15. li these evolutionary models can pro- 
duce magnetic fields strong enough to partially confine 
accreting nuclear fuel near the star's magnetic poles, even 
though the dipolc component of the magnetic field is very 
weak. This picture of magnetic field evolution in turn 
suggests that the shapes and phases of the nuclear- and 
accretion-powered pulses are similar to one another in 
some AMXPs because the nuclear- and accretion pow- 
ered X-ray emission comes from approximately the same 
area on the stellar surface. 

Effects of rapid spot motions. In Section 14.51 we 
discussed the effects of spot movements on timescales 
shorter than the time required to construct a pulse pro- 
file. Such effects are expected, because emitting spots 
are likely to move on the stellar surface on timescales 
at least as short as the ~0.1 ms dynamical time there 
whereas constructing a pulse profile usually requires fold- 
ing ~ 10 2 -10 3 s of X-ray flux data. Rapid spot mo- 
tions will produce X-ray flux variations on these same 
timescales. Our computations show that motion of the 



emitting area on the stellar surface on timescales longer 
than the spin period usually changes the amplitudes and 
the phases of the harmonic components of the theoretical 
pulse profile. 

Variations of the X-ray flux on any timescales shorter 
than the time required to construct a pulse profile will 
appear in the analysis as noise in excess of the normal 
counting noise. This noise will reduce the measured am- 
plitude of the oscillations at the spin frequency and its 
overtones for all spot locations, but its effect is likely to 
be stronger when the emitting area is near the spin axis 
because displacement of the emitting area by a given dis- 
tance there produces a larger change in the pulse phase 
and, for many geometries, in the pulse amplitude. This 
effect will therefore tend to reduce the apparent pulse 
amplitude even further when the emitting area is close 
to the spin axis. 

Undetectable and intermittent pulsations. In Sec- 
tion 13.11 we showed that if the emitting areas of some 
AMXPs are very close to the spin axis and remain there, 
the amplitudes of the oscillations they would produce can 
be ~0.5% or less, making them undetectable with cur- 
rent instruments. Rapid X-ray flux variations will make 
accretion-powered oscillations at the spin frequency more 
difficult to detect. Other effects, such as scattering of X- 
ray photons in circumstellar gas, may also play a role 
in reducing the detectability of such oscillations. These 
results show that the nearly aligned moving spot model 
may, possibly in combination with other effects, explain 
the nondetection of accretion-powered oscillations at the 
millisecond spin frequencies of some accreting neutron 
stars in which nuclear-powered oscillations have been de- 
tected. In Section 14.11 we showed that if the emitting 
area is within a few degrees of the spin axis and moves to- 
ward the rotation equator by ~ 10°, oscillations that were 
undetectable can become detectable. This may explain 
why accretion-powered oscillations appear only int ermit- 
tently in some AMXPs (|Lamb fc Bou tloukos 2008). The 
model suggests that oscillations may also disappear in- 
termittently in some AMXPs. 

Evolution of AMXP magnetic fields. In Section fe-ll we 
explained why the magnetic poles of most AMXPs are 
expected to be very close to their spin axes. One conse- 
quence is that their magnetic fields will channel accreting 
gas to the stellar surface near the spin axis. Hence the 
star's accretion-powered X-ray emission will come from 
areas near the spin poles. A second consequence is that 
many AMXPs may have surface magnetic fields as strong 
as ~ 10 11 -10 12 G, even though the dipole components 
of these fields are only ~ 10 8 -10 9 G. If so, many MRPs 
may also have surface magnetic fields as strong as ~ 10 11 - 
10 12 G, even though the dipole components inferred from 
their spin-down rates are only ~ 10 8 -10 9 G. 

Transient nature of AMXPs. In Section [5.21 we noted 
that the picture of AMXP magnetic field evolution just 
described suggests why the AMXPs in which accretion- 
powered oscillations have been detected arc in transient 
systems. If the magnetic poles of most neutron stars in 
LMXBs were forced very close to their spin axes during 
the initial, persistent phase of mass transfer, accretion- 
powered X-ray oscillations would be difficult or impos- 
sible to detect. Later, when mass transfer is transient, 
the stars will spin down, causing their magnetic poles to 
move away from their spin axes and making accretion- 
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powered X-ray oscillations detectable. 

Tests of the model. The nearly aligned moving spot 
model leads to a number of expectations about AMXP X- 
ray emission that can be tested (see Section[5]for details): 

1. The amplitudes, harmonic content, and arrival 
times of pulses should be functions of the X-ray 
luminosity and spectrum of the pulsar. 

2. The amplitudes, harmonic content, and arrival 
times of pulses should be correlated with one an- 
other. 

3. The arrival times of pulses with low amplitudes 
should fluctuate much more than the arrival times 
of pulses with high amplitudes. 

4. Pulse phase residuals are likely to form a track 
in the phase-residual versus pulse-amplitude plane. 
Such tracks are likely to be more evident if the 
range of the pulse amplitude variation is large. The 
position of pulses along such a track should be cor- 
related with the X-ray luminosity and spectrum of 
the pulsar. 

5. If the accretion- and nuclear-powered pulses of an 
AMXP appear nearly identical, long-term (days- 
to-weeks) variations in the phase residuals of the 
two types of pulses should track one another. 

6. The strength of the excess background noise pro- 
duced by pulse shape fluctuations should be corre- 
lated with the amplitudes, harmonic content, and 
arrival times of pulses and the X-ray luminosity 
and spectrum of the pulsar. 

7. The excess noise produced by motion of the emit- 
ting area should be stronger when the pulse ampli- 
tude is smaller and weaker when the pulse ampli- 
tude is larger. 

8. If AMXPs do have surface magnetic fields as 
strong as ~ 10 11 -10 12 G, their keV X-ray spec- 
tra may show strong- magnetic-field features. Such 
features are more likely in AMXPs that show 
evidence of nuclear fuel confinement, such as 
SAX J1808.4-3658 and XTE J1814-338. 

9. If AMXPs have surface magnetic fields ~ 10 11 - 
10 12 G, their offspring should have millisecond 
spin periods and total magnetic fields of simi- 
lar strength, but dipole fields - 10 8 -10 9 G. This 



should affect particle acceleration and 7-ray emis- 
sion by these neutron stars. 



10. If the explanation of the transient nature of the 
AMXPs suggested here is correct, they should be 
spinning down on long timescales. 



Note added in manuscript: This work was presented 
at the 2008 April Amsterdam Workshop on Accreting 
Millisecond Pulsars, where we emphasized a test of our 
nearly aligned moving spot model that could also have 
wider implications. We argued that the close similar- 
ity of the accretion-powered and burst oscillation wave- 
forms in XTE J1814— 338 strongly suggests that the 
emitting regions that produce them are similar and col- 
located, and that if, as we had proposed previously, the 
phase wandering of the accretion-powered oscillations is 
caused by movement of the emitting region, then the 
phase of the burst oscillatio ns should wander in the same 
way. After the Workshop, iWatts et ail (|2008l ) investi- 
gated this possibility and found just such a correlation 
in the RXTE data on XTE J1814-338. In particular, 
these authors found that during the main part of the two- 
month outburst of XTE J1814— 338, its burst oscillation 
not only has a wave form similar to that of the accretion- 
powered os cillation (|Strohmaver et al.1l2003t IWatts et al.l 
120051 [20061 ). but is also phase-locked with it, and that the 
peak of the burst oscillation coincides with the peak of 
the soft component of the accretion-powered oscillation. 
This indicates that, as we had suggested, the accretion- 
and nuclear-powered emitting regions in this pulsar very 
nearly coincide, and that the simultaneous wandering of 
the arrival times of both oscillations by ~ 1 ms (~ 0.3 in 
phase) during the outburst is due to wandering of the 
matter (and hence the fuel) deposition pattern on the 
stellar surface. 
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